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1.  INTRODUCTION 


Tuberous  sclerosis  complex  (TSC)  is  an  autosomal  dominant  disease  characterized  by  benign 
tumors  in  various  tissues.  The  genes  mutated  in  this  disease,  TSC1  and  TSC2,  encode  tumor 
suppressors  that  are  associated  in  a  complex.  The  TSC1/2  complex,  through  its  Rheb-GAP 
activity,  is  a  critical  negative  regulator  of  mTORCI  under  physiological  conditions.  Activation  of 
mTORCI  positively  stimulates  cap-dependent  mRNA  translation  via  its  downstream  substrates 
S6K  and  4E-BP.  In  our  previous  study,  we  demonstrated  that  TSC-mTORC1  signaling 
regulates  the  balance  between  cap-dependent  and  cap-independent  translation.  In  this  project, 
we  aim  to  (1)  Dissect  the  molecular  linkage  between  mTORCI  and  protein  homeostasis;  (2) 
Define  the  role  of  mTORCI  in  ribosome  dynamics  and  translational  re-programming. 
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2.  KEYWORDS 


Tuberous  sclerosis  complex;  mTORCI;  rapamycin;  mRNA;  translation;  ribosome;  quality 
control;  alternative  translation 
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3.  ACCOMPLISHMENTS 


•  What  were  the  major  goals  of  the  project? 

Task  1.  Dissect  the  Molecular  Linkage  between  mTORCI  and  Protein  Homeostasis 

Task  2.  Define  the  Role  of  mTORCI  in  Ribosome  Dynamics  and  Translational  Re¬ 
programming 

•  What  was  accomplished  under  these  goals? 

1.  TISdb:  a  database  for  alternative  translation  initiation  in  mammalian  cells. 

Proper  selection  of  the  translation  initiation  site  (TIS)  on  mRNAs  is  crucial  for  the  production  of 
desired  protein  products.  Recent  studies  using  ribosome  profiling  technology  uncovered  a 
surprising  variety  of  potential  TIS  sites  in  addition  to  the  annotated  start  codon.  The  prevailing 
alternative  translation  reshapes  the  landscape  of  proteome  in  terms  of  diversity  and  complexity. 
To  identify  the  hidden  coding  potential  of  transcriptome  in  mammalian  cells,  we  developed 
global  translation  initiation  sequencing  (GTI-Seq)  that  maps  genome-wide  TIS  positions  at 
nearly  single  nucleotide  resolution.  To  facilitate  studies  of  alternative  translation,  we  created  a 
database  of  alternative  TIS  sites  identified  from  human  and  mouse  cell  lines  based  on  multiple 
GTI-Seq  replicates.  The  TISdb,  available  at  http://tisdb.human.cornell.edu,  includes  6,991  TIS 
sites  from  4,961  human  genes  and  9,973  TIS  sites  from  5,668  mouse  genes.  The  TISdb 
website  provides  a  simple  browser  interface  for  query  of  high-confidence  TIS  sites  and  their 
associated  open  reading  frames  (ORFs).  The  output  of  search  results  provides  a  user-friendly 
visualization  of  TIS  information  in  the  context  of  transcript  isoforms.  Together,  the  information  in 
the  database  provides  an  easy  reference  for  alternative  translation  in  mammalian  cells  and  will 
support  future  investigation  of  novel  translational  products. 

2.  Ribosome  profiling  reveals  sequence-independent  post-initiation  pausing  as  a  signature  of 
translation. 


DMSO 


CHX 


The  journey  of  a  newly  synthesized  polypeptide  starts  in  the  peptidyltransferase  center  of  the 
ribosome,  from  where  it  traverses  the  exit  tunnel.  The  interior  of  the  ribosome  exit  tunnel  is 
neither  straight  nor  smooth.  How  the  ribosome  dynamics  in  vivo  is  influenced  by  the  exit  tunnel 
is  poorly  understood.  Genome-wide  ribosome  profiling  in  mammalian  cells  reveals  elevated 
ribosome  density  at  the  start  codon  and  surprisingly  the  downstream  5th  codon  position  as  well 
(Fig.  1 ).  We  found  that  the  highly  focused  ribosomal  pausing  shortly  after  initiation  is  attributed 
to  the  geometry  of  the  exit  tunnel,  as  deletion  of  the  loop  region  from  ribosome  protein  L4 
diminishes  translational  pausing 
at  the  5th  codon  position. 

Unexpectedly,  the  ribosome 
variant  undergoes  translational 
abandonment  shortly  after 
initiation,  suggesting  that  there 
exists  an  obligatory  step  between 
initiation  and  elongation 

commitment.  We  propose  that  pjg  ^etagene  analysis  of  RPFs  obtained  from  HEK293  cells  treated 

the  post-initiation  pausing  of  with  chx  (right)  or  dmso  control  (left). 
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ribosomes  represents  an  inherent  signature  of  the  translation  machinery  to  ensure  productive 
translation. 

3.  Quantitative  profiling  of  initiating  ribosomes  in  vivo. 

Cells  have  evolved  exquisite  mechanisms  to  fine-tune  the  rate  of 
protein  synthesis  in  response  to  stress.  Systemic  mapping  of 
start-codon  positions  and  precise  measurement  of  the 
corresponding  initiation  rate  would  transform  our  understanding 
of  translational  control.  Here  we  present  quantitative  translation 
initiation  sequencing  (QTI-seq),  with  which  the  initiating 
ribosomes  can  be  profiled  in  real  time  at  single-nucleotide 
resolution  (Fig.  2).  Resultant  initiation  maps  not  only  delineated 
variations  of  start-codon  selection  but  also  highlighted  a  dynamic 
range  of  initiation  rates  in  response  to  nutrient  starvation.  The 
integrated  data  set  provided  unique  insights  into  principles  of 
alternative  translation  and  mechanisms  controlling  different 
aspects  of  translation  initiation.  With  RiboTag  mice,  QTI-seq 
permitted  tissue-specific  profiling  of  initiating  ribosomes  in  vivo. 

Liver  cell-specific  ribosome  profiling  uncovered  a  robust 
translational  reprogramming  of  the  proteasome  system  in  fasted 
mice.  Qur  findings  illuminated  the  prevalence  and  dynamic  nature 
of  translational  regulation  pivotal  to  physiological  adaptation  in 
vivo. 

4.  Translational  control  of  cytosolic  stress  response  by  mitochondrial  ribosomal  protein  L18. 

In  response  to  stress,  cells  attenuate  global  protein  synthesis  but  permit  efficient  translation  of 
mRNAs  encoding  heat  shock  proteins  (HSPs).  Despite  decades  since  the  first  description  of  the 
heat  shock  response,  how  cells  achieve  translational  control  of  HSP  synthesis  remains 
enigmatic.  Here  we  report  an  unexpected  role  for  mitochondrial  ribosomal  protein  L18 
(MRPL18)  in  mammalian  cytosolic  stress  response.  MRPL18  bears  a  downstream  CUG  start 
codon  and  generates  a  cytosolic  isoform  in  a  stress-dependent  manner.  The  cytosolic  MRPL18 
incorporates  into  the  80S  ribosome  and  facilitates  ribosome  engagement  on  mRNAs  selected 
for  translation  during  stress.  MRPL18  knockdown  has  minimal  effects  on  mitochondria  function, 
but  substantially  dampens  cytosolic  HSP  expression  at  the  level  of  translation.  Qur  results 
uncover  a  hitherto  uncharacterized  stress  adaptation  mechanism  in  mammalian  cells,  which 
involves  formation  of  a  “hybrid”  ribosome  responsible  for  translational  regulation  of  cytosolic 

Fig.  3.  A  model  for  MRPL18  in 
translational  control  during 
stress  conditions.  In  stressed 
cells,  MRPL18  undergoes 
alternative  translation  from  the 
downstream  CUG  start  codon, 
which  results  in  a  cytoplasmic 
isoform  of  MRPL18  (red)  (2). 
MRPL18(cyto)  incorporates 
into  the  80S  ribosome 
complex,  facilitating  the 
engagement  of  mRNAs  highly 
expressed  under  stress,  such 
as  Hspala  (green)  (3). 


stress  response  (Fig.  3). 
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- 

Fig.  2.  Schematic  of  Ribo- 
seq  (left)  and  QTI-seq  (right) 
procedures. 
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•  What  opportunities  for  training  and  professional  development  has  the  project 
provided? 

Nothing  to  report 

•  How  were  the  results  disseminated  to  communities  of  interest? 

Nothing  to  report 

•  What  do  you  plan  to  do  during  the  next  reporting  period  to  accomplish  the  goals? 

1.  Eiucidating  the  roie  of  mRNA  modification  in  transiationai  controi 

N6-methyiadenosine  (m6A)  is  a  ubiquitous  modification  in  mRNA  and  other  RNAs  across  most 
eukaryotes.  For  many  years,  however,  the  exact  functions  of  m6A  were  mot  cieariy  understood. 
The  discovery  that  the  fat  mass  and  obesity-associated  protein  (FTO)  is  an  m6A  demethyiase 
indicates  that  this  modification  is  reversibie  and  dynamicaiiy  reguiated.  Recent  m6A  profiiing  in 
mouse  and  human  ceiis  further  showed  that  this  modification  is  a  widespread  mark  in  coding 
and  noncoding  RNA  transcripts,  it  is  iikeiy  that,  reversibie  RNA  methyiation,  anaiogous  to 
reversibie  DNA  and  histone  modifications,  may  affect  gene  expression,  in  coiiaborating  with 
Sarnie  Jeffery’s  iab  at  Corneii  Weiii  Medicai  Schooi,  we  are  investigating  the  roie  of  m6A 
modification  in  transiation  controi.  Our  preiiminary  resuits  indicate  that  stress-inducibie  Hsp70 
transiation  is  reguiated  by  m6A  modification.  By  coupiing  m6A  profiiing  with  QTi-seq,  we  aim  to 
eiucidate  the  mechanism  underiying  m6A  modification-mediated  transiationai  controi. 

2.  Linking  mitochondriai  tRNA  modification  and  the  proteome  iandscape 

Mitochondrion  represents  an  eiegant  exampie  of  a  seif-sustained  system.  Bearing  its  own 
genome  and  the  transiation  machinery,  mitochondria  produces  proteins  necessary  for  ATP 
production.  Approximateiy  200  pathogenic  mutations  have  been  mapped  to  mitochondriai  tRNA 
(mt  tRNA)  genes.  Chemicai  modifications  are  a  characteristic  structurai  feature  of  mt  tRNAs 
which  contain  16  species  of  modified  nucieosides  inciuding  three  that  are  mitochondria-specific. 
Loss  of  modifications  found  at  the  wobbie  positions  of  mt  tRNAs  have  been  impiicated  in 
disrupted  tRNA  structures  that  may  affect  their  transiation  capabiiity.  However  it  is  stiii  not 
compieteiy  understood  how  mt  tRNA  modification  contributes  to  the  reguiation  of  the 
mitochondriai  proteome.  in  particuiar,  there  is  iittie  mechanistic  knowiedge  about  the  nature  of 
crosstaik  between  mitochondriai  activity  and  ceiiuiar  stress  signais. 

To  address  this  compiex  bioiogicai  probiem  in  a  systematic  manner,  we  deveioped  innovative 
approaches  by  coupiing  high-resoiution  ribosome  profiiing  with  quantitative  mapping  of  tRNA 
modification.  Mridu  Saikia,  the  new  postdoc  in  my  iaboratory,  hypothesized  that  mt  tRNA 
modification  acts  as  a  sentinei  for  ceiiuiar  stress  signais  and  the  aitered  proteome  triggers 
mitonuciear  protein  imbaiance  and  contributes  to  disease  deveiopment.  We  wiii  determine  how 
ceiiuiar  stress  affects  mitochondriai  metaboiism,  map  mt  tRNA  modification  in  response  to 
stress,  and  define  the  roie  of  mt  tRNA  modification  in  transiation. 
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4.  IMPACT 


•  What  was  the  impact  on  the  development  of  the  principal  discipline(s)  of  the  project? 

To  maintain  protein  homeostasis,  cells  have  evolved  several  highly  coordinated  systems  to 
sense  and  respond  to  environmental  cues.  Coupling  with  the  chaperone  network,  autophagy, 
and  ubiquitin/proteasome  system,  translation  lies  at  the  heart  of  protein  homeostasis. 
Translational  re-programming  allows  selective  translation  of  mRNAs  to  maintain  the  expression 
of  specific  proteins  during  patho-physiological  conditions.  Such  regulation  can  be  quantitative 
(all-or-none  vs.  graded),  qualitative  (enabling  a  single  mRNA  to  produce  several  different 
proteins),  or  selective  (activating  subsets  of  mRNAs  for  translation).  We  developed  quantitative 
profiling  of  initiating  ribosomes  that  permits  high  resolution  mapping  of  alternative  translation 
initiation  sites.  Our  study  reveals  that  translational  re-programming  plays  a  crucial  role  in  cellular 
metabolism. 

•  What  was  the  impact  on  other  disciplines? 

Nothing  to  report 

•  What  was  the  impact  on  technology  transfer? 

Nothing  to  report 

•  What  was  the  impact  on  society  beyond  science  and  technology? 

Nothing  to  report 
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5.  CHANGES  /  PROBLEMS 

Nothing  to  report 
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6.  PRODUCTS 


•  Publications,  conference  papers,  and  presentations 

Journal  publications  (*  Corresponding  author) 

1.  Gao  X,  Wan  J,  and  Qian  SB*.  Genome-wide  profiling  of  alternative  translation  initiation  sites. 
Methods  Mol  Biol  2015  (in  press) 

2.  Zhang  X,  Gao  X,  Roots  RA,  Conn  CS,  Liu  B,  and  Qian  SB*.  Translational  control  of  cytosolic 
stress  proteins  by  mitochondrial  ribosomal  protein  LI  8.  Nat  Struct  Mol  Biol  2015; 
22(5):404-10 

3.  Gao  X,  Wan  J,  Liu  B,  Ma  M,  Shen  B,  and  Qian  SB*.  Quantitative  profiling  of  initiating 
ribosomes  in  vivo.  Nat  Methods  2015;  12(2):  147-53 

4.  Han  Y,  Gao  X,  Liu  B,  Wan  J,  Zhang  X,  and  Qian  SB*.  Ribosome  profiling  reveals  sequence- 
independent  post-initiation  pausing  as  a  signature  of  translation.  Cell  Res  2014;  24(7):842- 
51 
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Acknowledgement  of  federal  support:  Yes. 


•  Other  products 

Databases:  http://tisdb.human.cornell.edu 

TISdb:  a  database  for  alternative  translation  initiation  in  mammalian  cells. 
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7.  PARTICIPANTS  &  OTHER  COLLABORATING  ORGANIZATIONS 


•  What  individuals  have  worked  on  the  project? 

Name:  Shu-Bing  Qian,  PhD 

Project  Role:  PI 

Nearest  person  month  worked:  1 

Contribution  to  Project:  Dr.  Qian  is  responsible  for  the  overall  administration  and 

direction  of  the  project,  including  designing  experimental  protocols,  interpreting  results,  writing 
manuscripts  and  supervising  research  assistants. 


Name: 

Project  Role: 

Nearest  person  month  worked: 
Contribution  to  Project: 
experiments  (Aim  1,  Aim  2). 


Jun  Zhou,  PhD 
Postdoc 
12 

Dr.  Zhou  will  be  responsible  for  ribosome  profiling 


•  Has  there  been  a  change  in  the  active  other  support  of  the  PD/PI(s)  or  senior/key 
personnel  since  the  last  reporting  period? 

Nothing  to  report 

•  What  other  organizations  were  involved  as  partners? 

Nothing  to  report 
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©2014  IBCB.SIBS,  CAS  All  rights  reserved  1001-0602/14  $  32.00  ' 

ORIGINAL  ARTICLE  www.nature.com/cr 

Ribosome  profiling  reveals  sequence-independent 
post-initiation  pausing  as  a  signature  of  translation 

Yan  Xiangwei  Gao'  *,  Botao  Liu\  Ji  Wan',  Xingqian  Zhang',  Shu-Bing  Qian'’^ 

'Division  of  Nutritional  Sciences,  Cornell  University,  Ithaca,  NY  14853,  USA;  ^Department  of  Infectious  Diseases,  Ruijin  Hospital, 
Shanghai  Jiaotong  University  School  of  Medicine,  Shanghai  200025,  China;  ^Graduate  Field  of  Genetics,  Genomics  and  Develop¬ 
ment,  Cornell  University,  Ithaca,  NY  14853,  USA 


The  journey  of  a  newly  synthesized  polypeptide  starts  in  the  peptidyltransferase  center  of  the  rihosome,  from 
where  it  traverses  the  exit  tunnel.  The  interior  of  the  rihosome  exit  tunnel  is  neither  straight  nor  smooth.  How  the 
rihosome  dynamics  in  vivo  is  influenced  hy  the  exit  tunnel  is  poorly  understood.  Genome-wide  rihosome  profiling  in 
mammalian  cells  reveals  elevated  rihosome  density  at  the  start  codon  and  surprisingly  the  downstream  5th  codon  po¬ 
sition  as  well.  We  found  that  the  highly  focused  rihosomal  pausing  shortly  after  initiation  is  attributed  to  the  geome¬ 
try  of  the  exit  tunnel,  as  deletion  of  the  loop  region  from  rihosome  protein  L4  diminishes  translational  pausing  at  the 
5th  codon  position.  Unexpectedly,  the  rihosome  variant  undergoes  translational  abandonment  shortly  after  initiation, 
suggesting  that  there  exists  an  obligatory  step  between  initiation  and  elongation  commitment.  We  propose  that  the 
post-initiation  pausing  of  ribosomes  represents  an  inherent  signature  of  the  translation  machinery  to  ensure  produc¬ 
tive  translation. 

Keywords',  ribosome;  translation 

Cell  Research  (2014)  24:842-851.  doi:  10. 1038/cr.2014.74;  published  online  6  June  2014 


Introduction 

Ribosomes  decode  the  genetic  information  in  mRNA 
and  convert  it  into  the  amino  acid  sequences  of  proteins. 
During  mRNA  translation,  the  ribosome  does  not  pro¬ 
ceed  at  a  constant  rate  but  rather  in  a  stop-and-go  move¬ 
ment  manner  [1].  Variations  of  elongation  speed  may 
result  from  local  stable  mRNA  structures  [2],  the  pres¬ 
ence  of  rare  codons  [3],  or  the  interactions  between  the 
nascent  chain  and  the  ribosome  and/or  exogenous  factors 
[4].  Frequent  translational  pausing  may  have  functional 
importance  for  the  cellular  production  line.  For  instance, 
ribosome  stalling  could  trigger  frame-shifting  or  transla¬ 
tional  abandonment  [5,  6].  A  growing  body  of  evidence 
suggests  that  cells  purposely  employ  ribosomal  pausing 
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as  a  regulatory  mechanism  to  maintain  intracellular  pro¬ 
tein  homeostasis.  For  example,  translational  delays  in 
specific  regions  facilitate  co-translational  folding  of  na¬ 
scent  chains  [7].  Ribosomes  also  adopt  early  elongation 
pausing  in  response  to  various  stress  conditions  [8,  9]. 
Given  the  wide  range  of  impacts  of  ribosome  pausing  on 
cellular  systems,  a  fundamental  question  is  to  dissect  the 
mechanistic  basis  underlying  translational  pausing  and 
determine  the  physiological  consequences  when  ribo¬ 
somes  stall  on  individual  mRNAs. 

Despite  several  decades  since  the  first  description  of 
ribosome  pausing  and  stacking  during  translation  [10], 
it  was  not  possible  until  now  to  measure  ribosome  dy¬ 
namics  in  vivo  at  the  genome-wide  scale.  Ribosome 
profiling  or  Ribo-seq,  based  on  deep  sequencing  of  ribo¬ 
some-protected  mRNA  fragments  (RPFs),  has  proven  to 
be  powerful  in  determining  ribosome  positions  and  den¬ 
sities  across  the  entire  transcriptome  [11].  Many  research 
groups  have  successfully  applied  Ribo-seq  to  a  wide 
range  of  organismal  species  and  reported  pervasive  pause 
sites  along  the  coding  region  (CDS)  of  mRNAs  [12-18]. 
From  bacteria  to  human  cells,  the  most  consistent  pause 
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sites  are  located  at  the  start  and  stop  codons,  which  is  in 
line  with  the  notion  that  ribosomes  tend  to  pause  during 
initiation  and  termination  [10].  Ribosome  profiling  also 
captures  many  sequence-dependent  stalling  sites  [15]. 
However,  much  less  is  known  about  translational  pausing 
governed  by  factors  other  than  cw-elements. 

Here,  we  report  and  characterize  a  specific  ribosomal 
pausing  at  the  5th  codon  position  of  mRNA  CDSs  in 
mammalian  cells.  Unlike  sequence-dependent  elongation 
pausing,  the  ribosome  stalling  shortly  after  initiation  is 
attributed  to  the  geometry  of  the  exit  tunnel.  Importantly, 
modifying  the  exit  tunnel  not  only  reduces  the  ribosomal 
pausing  at  the  5th  codon  position,  but  results  in  an  unex¬ 
pected  translational  abandonment.  Our  finding  suggests 
a  rich  communication  between  the  nascent  chain  and  the 
translation  machinery  in  regulating  ribosome  kinetics 


during  translation. 

Results 

Ribosome  profiling  reveals  post-initiation  pausing 

We  applied  a  modified  Ribo-seq  to  HEK293  cells  and 
examined  the  average  pattern  of  ribosome  density  across 
the  entire  transcriptome  [8,  19].  As  expected,  the  mapped 
RPF  reads  start  an  abrupt  accumulation  at  the  annotated 
start  codon  [11]  (Figure  lA).  Interestingly,  we  observed 
another  prominent  peak  at  the  12  nt  downstream  posi¬ 
tion,  corresponding  to  the  5th  codon  position  of  the  CDS 
(Figure  lA,  red  arrow).  This  is  not  due  to  cycloheximide 
(CHX)-mediated  translational  arrest  because  the  same 
feature  is  maintained  in  cells  without  any  dmg  treatment 
(Figure  lA,  left  panel).  The  ribosomal  pausing  at  the 
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Figure  1  Characterization  of  post-initiation  pausing.  (A)  Metagene  analysis  of  RPFs  obtained  from  HEK293  cells  treated  with 
CHX  (right)  or  DMSO  control  (left).  All  mapped  reads  are  aligned  at  the  annotated  start  codon  AUG  and  stratified  by  frames. 
The  read  density  at  each  nucleotide  position  is  averaged  using  the  P-site  of  RPFs.  The  red  arrow  indicates  the  5th  codon 
position.  (B)  Read  fractions  relative  to  the  annotated  reading  frame  are  calculated  for  RPFs  mapped  to  the  1st  (AUG)  codon, 
5th  codon,  or  the  entire  CDS.  Data  are  shown  as  means  ±  SD,  n  =  4.  (C)  Metagene  analysis  of  RPFs  obtained  from  HEK293 
cells  treated  with  LTM  (upper)  or  CHX  (bottom).  All  mapped  reads  are  aligned  at  the  annotated  start  codon  (aTIS,  green  tri¬ 
angle;  left  panel),  upstream  TIS  (uTIS,  white  triangle;  middle  panel),  or  downstream  TIS  (dTIS,  blue  triangle;  right  panel).  The 
read  density  at  each  nucleotide  position  is  averaged  using  the  P-site  of  RPFs.  The  red  arrow  indicates  the  5th  codon  position. 


www.cell-research.com  |  Cell  Research 


Ribosomal  pausing  before  elongation  commitment 

m 


5th  codon  position  is  highly  reproducible  as  evidenced 
by  the  similar  pattern  of  RPF  distribution  in  a  total  of  4 
independent  biological  replicates  (Supplementary  infor¬ 
mation,  Figure  SI).  To  confirm  whether  the  same  feature 
also  exists  in  different  cell  types,  we  conducted  Ribo-seq 
in  FleLa  cells  and  a  mouse  embryonic  fibroblast  (MEF) 
cell  line.  Both  cells  showed  consistent  pattern  of  RPFs 
with  prominent  peaks  of  ribosome  density  at  both  the 
start  codon  and  the  5th  codon  position  (Supplementary 
information.  Figure  S2).  Thus,  the  excess  of  ribosome 
density  at  the  5th  codon  position  represents  a  strong  ribo¬ 
somal  pausing  shortly  after  initiation. 

The  post-initiation  pausing  bears  several  unique  fea¬ 
tures  that  distinguish  it  from  the  typical  ribosome  paus¬ 
ing  at  the  start  codon.  First,  the  RPFs  recovered  at  the 
5th  codon  position  demonstrate  a  near  perfect  framing 
or  phasing  with  >  90%  of  reads  mapped  to  the  same 
frame  as  the  annotated  ORFs  (Figure  lA  and  IB).  This 
is  in  contrast  to  the  RPF  reads  derived  from  other  CDS 
regions  including  the  AUG  start  codon,  in  which  about 
a  quarter  of  the  reads  are  out  of  the  annotated  reading 
frame  (Figure  IB).  Second,  RPFs  derived  from  the  5th 
codon  position  have  different  size  distribution  in  length 
(Supplementary  information.  Figure  S3).  While  the  ma¬ 
jority  of  footprints  mapped  to  the  entire  CDS  were  29- 
mer,  a  substantial  proportion  of  RPFs  mapped  to  the  5th 
codon  position  are  27-mer.  These  results  suggest  that  the 
stalled  ribosome  at  the  5th  codon  position  might  adopt  a 
more  compact  conformation  than  the  others.  Supporting 
this  notion,  recent  studies  demonstrated  that  the  length  of 
footprints  could  be  infiuenced  by  ribosome  dynamics  [20, 
21]. 

To  rule  out  the  possibility  that  the  prominent  ribosom¬ 
al  pausing  at  the  5th  codon  position  is  due  to  reinitiation, 
we  segregated  initiating  ribosomes  from  elongating 
ribosomes  by  applying  a  translation  inhibitor  lactimido- 
mycin  (LTM).  Unlike  CFIX  that  freezes  all  the  ribosomes 
engaged  on  mRNAs,  LTM  preferentially  acts  on  the 
initiating  ribosome  but  not  the  elongating  one  [22].  As 
we  reported  previously,  LTM-based  ribosome  profiling 
(GTI-seq)  offers  global  mapping  of  translation  initiation 
sites  (TIS)  at  single-nucleotide  resolution  [19].  GTI- 
seq  revealed  a  pronounced  single  peak  highly  focused 
at  the  annotated  TIS  wilh  few  reads  recovered  in  Ihe 
downstream  region  (Figure  1C,  left  panel).  The  lack  of 
LTM-associated  RPF  reads  at  the  5th  codon  position 
strongly  indicates  that  the  ribosomal  pausing  at  this  po¬ 
sition  is  a  true  elongation  event.  GTI-seq  also  permits 
identification  of  alternative  start  codons,  including  up¬ 
stream  TIS  and  downstream  TIS  sites  [19].  To  examine 
whether  alternative  initiation  is  also  followed  by  down¬ 
stream  ribosomal  pausing,  we  aligned  transcripts  by  the 


alternative  TIS  sites  idenlified  by  GTI-seq  (Figure  1C, 
middle  and  righi  panels).  Remarkably,  the  CFIX-associ- 
ated  RPFs  exhibit  dominant  peaks  not  only  at  the  corre¬ 
sponding  TIS  posilions,  but  also  at  the  downstream  5th 
codon  position.  Thus,  post-initiation  pausing  represents  a 
universal  feature  for  ribosomes  after  initiation  from  both 
the  annotated  and  the  alternative  start  codons. 

When  the  ribosomal  P-site  is  mapped  at  the  5th  codon 
position,  its  footprint  coincidently  starts  with  AUG  at  the 
5'  end.  One  possibility  is  that  these  RPFs  are  preferen¬ 
tially  obtained  after  RNase  I  digestion  or  selectively  am¬ 
plified  during  deep  sequencing.  Flowever,  this  is  not  the 
case.  The  frequency  of  RPFs  with  AUG  at  the  5'  end  is 
among  the  average  compared  to  RPFs  starting  with  other 
sequences  (Supplementary  information.  Figure  S4).  Fur¬ 
thermore,  many  internal  AUG  codons  are  not  associated 
with  the  downstream  pausing.  Therefore,  the  elevated 
ribosome  density  at  the  5th  codon  position  after  the  an¬ 
notated  start  codon  is  likely  coupled  with  the  initiation 
event. 

Characterization  of  post-initiation  pausing 

To  investigate  possible  mechanisms  underlying 
post-initiation  pausing  at  the  5th  codon  position,  we  first 
examined  whether  ribosome  stalling  at  this  specific  posi¬ 
tion  is  sequence  dependent.  Flowever,  we  found  a  similar 
codon  composition  at  the  5th  codon  position  to  the  one 
averaged  across  the  entire  CDS  (Supplementary  informa¬ 
tion,  Figure  S5).  To  find  possible  general  relationships 
between  the  functions  of  genes  and  post-initiation  paus¬ 
ing,  we  searched  for  common  biological  themes  in  tran¬ 
scripts  with  high  and  low  ribosomal  pausing  at  the  5th 
codon  position.  Gene  ontology  term  analysis  revealed  no 
specific  gene  functions  associated  with  either  strong  or 
weak  post-initiation  pausing  (Supplementary  informa¬ 
tion,  Figure  S6).  These  results  indicate  that  the  ribosomal 
pausing  shortly  after  initiation  is  sequence  independent, 
likely  representing  a  translational  event  inherently  asso¬ 
ciated  with  the  translation  machinery  per  se. 

We  next  asked  whether  ribosomal  pausing  at  the  5th 
codon  position  is  a  result  of  ribosomal  stacking  during 
elongation.  To  this  end,  we  isolated  the  monosome  frac¬ 
tion  followed  by  Ribo-seq  (Supplementary  information. 
Figure  S7).  Since  every  transcript  in  the  monosome  has 
only  one  ribosome,  the  ribosome  behavior  is  not  infiu¬ 
enced  by  any  neighboring  ribosomes.  Despite  much  less 
RPF  density  in  the  CDS  region,  the  monosome  fraction 
maintained  prominent  peaks  at  both  the  start  codon  and 
the  5th  codon  position  (Supplementary  information. 
Figure  S7).  This  result  further  supports  the  notion  that 
post-initiation  pausing  reflects  an  intrinsic  feature  of  ri¬ 
bosomes  shortly  after  initiation. 
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Having  characterized  the  post-initiation  pausing  in 
mammalian  cells,  we  attempted  to  expand  our  findings 
to  other  species  whose  Ribo-seq  data  are  available  in  the 
literature.  Reanalyzing  the  RPF  reads  using  the  same 
pipeline  revealed  similar  post-initiation  pausing  in  other 
eukaryotic  cells  (Supplementary  information,  Figure  S8), 
including  budding  yeast  S.  cerevisiae  [11],  nematode  C. 
elegans  [17],  and  zebrafish /).  rerio  [16].  Interestingly, 
the  most  recent  report  of  mitochondrial  ribosome  profil¬ 
ing  also  showed  similar  ribosomal  pausing  at  the  5th  co¬ 
don  position  [23].  Notably,  the  relative  extent  of  post-ini¬ 
tiation  pausing  varies  among  different  studies  and  we 
cannot  exclude  the  possibility  that  technical  differences 
(e.g.,  Mg^^  concentration  in  polysome  buffer)  may  cause 
variations  of  ribosome  density  at  the  5th  codon  position. 
Considering  its  omnipresence  in  multiple  species,  it  is 
likely  that  the  post-initiation  pausing  is  a  general  feature 
for  different  types  of  translation  machinery. 

Post-initiation  pausing  is  attributed  to  the  geometry  of 
the  exit  tunnel 

Given  the  sequence-independent  but  position-specif¬ 
ic  feature,  the  observed  post-initiation  pausing  is  likely 
governed  by  the  ribosome  exit  tunnel.  The  interior  of  the 
ribosome  exit  tunnel  is  neither  straight  nor  smooth  with 
an  evident  constriction  site  formed  by  RPL4  and  RPL17 
[24]  (Figure  2A).  To  investigate  whether  the  highly 
focused  ribosomal  pausing  shortly  after  initiation  is  at¬ 
tributed  to  the  geometry  of  the  exit  tunnel,  we  focused 
on  RPL4  because  it  is  the  first  ribosomal  protein  that  the 
nascent  chain  potentially  interacts  within  the  exit  tunnel. 
We  constructed  a  RPL4  mutant  by  deleting  the  most  con¬ 
served  HRSG  loop  region  [25]  (Supplementary  informa¬ 
tion,  Figure  S9).  After  overexpression  in  HEK293  cells, 
myc-tagged  RPL4(Aloop)  mutant  showed  a  similar  cellu¬ 
lar  localization  to  the  wild  type  (Figure  2B).  In  addition 
to  the  nucleolus  localization  that  is  typical  for  ribosome 
subunit  assembly,  both  the  wild-type  and  the  mutant 
RPL4  showed  the  same  cytoplasmic  distribution  as  ribo¬ 
somal  PO,  a  protein  associated  only  with  the  assembled 
ribosome  complex.  To  confirm  that  RPL4(Aloop)  is  able 
to  incorporate  into  the  assembled  SOS  ribosome  complex 
in  the  cytoplasm,  we  performed  immunoprecipitation 
(IP)  assay  using  an  anti-myc  antibody  in  HEK293  cells 
stably  expressing  wild-type  RPL4  or  RPL4(Aloop).  The 
presence  of  a  myc  tag  at  the  C  terminus  of  the  transgene 
prevents  precipitation  of  ribosomes  synthesizing  the 
RPL4  polypeptide  itself.  In  addition,  we  converted  the 
polysome  into  monosome  by  RNase  I  digestion  prior  to 
anti-myc  IP  to  exclude  indirect  ribosome  pull-down  (Fig¬ 
ure  2C).  Both  small  and  large  ribosomal  proteins  were 
readily  precipitated  by  anti-myc  antibody  from  cells  ex- 
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Figure  2  Reduced  post-initiation  pausing  for  ribosomes  bearing 
L4  mutant.  (A)  Schematic  structure  of  ribosomes  highlighting 
the  exit  tunnel.  tRNA  and  the  nascent  chain  are  shown  in  or¬ 
ange,  RPL4  in  red,  and  RPL17  in  blue.  (B)  Cellular  localization 
of  L4  wild  type  and  Aloop  mutant  is  examined  in  HEK293  cells 
using  immunostaining.  Scale  bar,  10  pm.  (C)  Myc-tagged  L4 
is  incorporated  into  ribosome  complexes.  Polysome  fractions 
from  transfected  HEK293  cells  were  treated  with  RNase  I  prior 
to  anti-myc  IP.  Both  the  input  and  the  immunoprecipitates  were 
blotted  using  antibodies  as  indicated.  (D)  Metagene  analysis 
of  RPFs  obtained  from  purified  ribosomes  bearing  L4  wild  type 
(top)  and  Aloop  mutant  (bottom).  All  mapped  reads  are  aligned 
at  the  annotated  start  codon  AUG  and  stratified  by  frames.  (E) 
Read  density  over  CDS  relative  to  the  initiation  site  (1st  -  5th)  is 
shown  in  box  plots. 


pressing  either  wild-type  RPL4  or  RPL4(Aloop)  (Figure 
2C,  lanes  5  and  6),  but  not  GFP  control.  Thus,  deleting 
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the  loop  region  does  not  prevent  RPL4  from  integrating 
into  the  assembled  SOS  ribosome  complexes. 

We  next  examined  the  behavior  of  ribosomes  bearing 
myc-tagged  RPL4  by  Ribo-seq  of  transfected  cells.  To 
exclude  the  RPFs  derived  from  ribosomes  containing 
endogenous  RPL4,  we  conducted  selective  ribosome 
profiling  by  purifying  myc-tagged  ribosomes  after  RNase 
I  digestion.  Thus,  all  the  obtained  RPFs  are  exclusively 
derived  from  ribosomes  bearing  myc-tagged  RPL4.  As 
expected,  ribosomes  containing  RPL4(WT)  showed  ele¬ 
vated  ribosome  density  at  both  the  start  codon  and  the  5th 
codon  positions  (Figure  2D).  Flowever,  RPL4(Aloop)-in- 
corporated  ribosomes  exhibited  a  substantial  reduction 
of  ribosome  density  at  the  5th  codon  position  relative  to 
the  start  codon  (Figure  2D,  bottom  panel).  The  dimin¬ 
ished  post-initiation  pausing  is  also  evident  on  individual 
transcripts,  including  the  relatively  long  P4HB  and  the 
short  RPL2 3  (Supplementary  information.  Figure  SIO). 
Notably,  the  ribosome  density  at  the  start  codon  remains 
high  for  ribosomes  containing  RPL4(Aloop),  strongly 
indicating  that  the  assembly  of  SOS  ribosomes  at  the  start 
codon  is  not  affected  by  this  ribosome  variant. 

Ribosomes  with  altered  exit  tunnel  lead  to  abortive  trans¬ 
lation 

For  ribosomes  with  modified  exit  tunnel,  the  reduced 
post-initiation  pausing  is  expected  to  promote  the  rate 
of  translation.  Yet  the  RPL4(Aloop)-associated  RPFs 
showed  an  obvious  decrease  of  footprint  density  over 
the  downstream  CDS  region  relative  to  the  start  codon 
(Figure  2E).  This  surprising  finding  suggests  that  the 
ribosome  variant  with  a  modified  exit  tunnel  does  not 
support  continuous  translation  elongation.  This  could 
be  explained  by  ribosome  drop-off  that  occurs  in  the 
presence  of  RPL4(Aloop).  To  test  this  possibility,  we  in¬ 
spected  the  distribution  of  myc-tagged  RPL4  in  ribosome 
fractions  resolved  by  sucrose  gradient  velocity  sedimen¬ 
tation  (Figure  3A).  Myc-tagged  RPL4(WT)  completely 
behaved  like  the  endogenous  one  as  evidenced  by  their 
nearly  identical  distribution  in  polysome  fractions. 
RPL4(Aloop),  in  contrast,  showed  an  evident  decrease  in 
the  polysome  fractions  relative  to  the  endogenous  RPL4 
despite  its  comparable  enrichment  in  the  monosome 
(Figure  3A  and  3B).  These  results  suggest  that  an  altered 
ribosomal  exit  tunnel  may  cause  ribosome  dissociation 
shortly  after  initiation. 

Notably,  the  distance  between  the  peptidyltransferase 
center  (PTC)  and  the  RPL4  loop  region  is  more  than  five 
amino  acids  [24].  Removing  the  RPL4  loop  region  does 
not  seem  to  clog  or  close  the  exit  tunnel  as  we  observed 
no  additional  pausing  sites  within  the  30-codon  window 
for  ribosomes  bearing  RPL4(Aloop)  (Figure  2D).  To  ex¬ 


amine  whether  translational  abandonment  mainly  occurs 
before  the  nascent  chain  emerges  out  of  the  ribosome,  we 
performed  nascent  chain  IP  in  HEK293  cells  expressing 
a  GFP  reporter  with  an  N-terminal  Flag-tag  [8].  Early 
ribosome  dissociation  would  prevent  anti-Flag  antibody 
from  precipitating  translating  ribosomes.  For  cells  ex¬ 
pressing  RPL4(WT),  anti-Flag  IP  could  pull  down  con¬ 
siderable  amount  of  myc-tagged  RPL4  along  with  other 
ribosome  proteins  (Figure  3C,  lane  5).  Flowever,  in  cells 
expressing  RPL4(Aloop),  very  little  RPL4(Aloop)  was 
co-precipitated  by  anti-Flag  antibody  when  compared  to 
other  ribosome  proteins  (Figure  3C,  lane  6).  This  result 
is  consistent  with  early  abandonment  of  translation  for 
ribosomes  with  modified  exit  tunnel. 

For  ribosomes  with  altered  exit  tunnel,  it  is  unclear 
which  event  happens  first:  reduced  post-initiation  paus¬ 
ing  or  ribosome  dissociation.  Although  the  lack  of 
post-initiation  pausing  could  trigger  ribosome  drop-off, 
it  is  equally  possible  that  the  diminished  post-initiation 
pausing  is  a  direct  consequence  of  abolished  translation 
shortly  after  initiation.  In  the  latter  scenario,  dissociation 
of  ribosomes  bearing  RPL4(Aloop)  must  occur  before 
the  nascent  chain  reaches  the  5th  codon  position.  Despite 
lots  of  efforts,  we  were  unable  to  determine  the  length  of 
tRNA-linked  short  peptides  associated  with  ribosomes 
bearing  RPL4(Aloop).  Nevertheless,  our  results  argue 
that  the  very  first  few  amino  acids  are  able  to  sense  the 
integrity  of  the  ribosomal  exit  tunnel  in  an  active  manner 
and  post-initiation  pausing  of  ribosomes  may  facilitate 
elongation  commitment. 

Ribosomes  with  altered  exit  tunnel  lead  to  cellular  toxicity 

The  frequent  translational  abandonment  is  expected  to 
be  detrimental  to  cellular  functions.  This  was  indeed  the 
case.  We  examined  cell  viability  after  knocking  down  the 
endogenous  RPL4  using  lentiviruses  expressing  small 
hairpin  RNA  (shRNA)  targeting  the  3'  UTR  of  RPL4. 
Consistent  with  the  essential  role  of  RPL4  in  protein 
production,  we  observed  progressive  cellular  toxicity 
in  control  cells  expressing  GFP  after  RPL4  knockdown 
(Figure  4A).  Only  cells  expressing  wild-type  RPL4,  but 
not  RPL4(loop),  were  able  to  maintain  higher  cell  via¬ 
bility.  Thus,  the  post-initiation  pausing  governed  by  the 
ribosomal  exit  tunnel  represents  an  obligatory  step  be¬ 
tween  initiation  and  elongation. 

Discussion 

A  longstanding  question  in  translational  control  is  to 
systemically  understand  physiological  purposes  of  ribo¬ 
somal  pausing.  mRNA  structure,  tRNA  abundance,  and 
amino  acid  sequence  all  have  the  potential  to  impact  the 
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Figure  3  Ribosomes  bearing  L4  mutant  lead  to  abortive  translation.  (A)  Polysome  profiles  of  HEK293  cells  expressing  L4 
wild  type  or  loop  mutant.  The  bottom  panel  shows  the  distribution  of  exogenous  and  endogenous  ribosome  proteins  in  ribo¬ 
some  fractions.  (B)  A  simple  model  of  translational  abandonment  for  ribosomes  bearing  L4(Aloop)  mutant.  The  right  panel 
shows  the  quantification  of  myc-tagged  L4  relative  to  the  endogenous  L4  in  ribosome  fractions  as  shown  in  A.  Data  are 
shown  as  means  ±  SD,  n  =  3.  (C)  Schematic  for  nascent  chain  IP  assay  using  a  HEK293  cell  line  expressing  flag-GFP  (top 
panel).  Transfected  HEK293  cells  were  lysed  and  treated  with  RNase  I  prior  to  anti-flag  IP.  Both  the  input  and  the  immuno- 
precipitates  were  blotted  for  ribosomal  proteins  as  indicated  (bottom  panel). 


ribosome  dynamics.  Based  on  multiple  empirical  data 
sets,  here  we  report  a  sequence-independent  but  posi¬ 
tion-specific  ribosomal  pausing  governed  by  the  geom¬ 
etry  of  the  exit  tunnel.  Several  lines  of  evidence  support 
the  notion  that  the  ribosome  pausing  at  the  5th  codon 
position  serves  as  a  functional  signature  of  translation. 
First,  post-initiation  pausing  of  ribosomes  is  nearly  a  uni¬ 
versal  feature  in  a  wide  range  of  organisms.  With  an  ever 
growing  application  of  ribosome  profiling  in  different 
cell  types  and  tissues,  we  expect  to  see  a  more  compre¬ 
hensive  picture  illuminating  ribosome  dynamics  between 
initiation  and  elongation.  Second,  we  provide  direct  ev¬ 
idence  that  modifying  the  exit  tunnel  infiuences  the  ex¬ 
tent  of  ribosomal  pausing  shortly  after  initiation.  To  our 
knowledge,  this  is  the  first  example  of  ribosomal  pausing 
associated  with  the  translation  machinery.  Third,  post-ini¬ 


tiation  pausing  of  ribosomes  is  of  physiological  signifi¬ 
cance  as  elimination  of  post-initiation  pausing  is  coupled 
with  translational  abandonment.  Together,  these  results 
argue  that  the  ribosome  behavior  is  intimately  controlled 
via  the  constant  interaction  between  the  growing  nascent 
chain  and  the  translation  machinery  (Figure  4B). 

It  has  long  been  believed  that,  after  the  assembly  of 
SOS  ribosome  at  the  start  codon,  elongation  ensues  via 
a  simple  repetitive  process  comprising  ribosomal  trans¬ 
location  and  tRNA  recycling.  However,  a  growing  body 
of  evidence  suggests  that  the  synthesis  of  the  very  first  4 
amino  acids  possesses  special  features.  In  prokaryotes, 
the  presence  of  macrolides  allows  synthesis  of  at  least  4 
amino  acids  before  peptidyl-tRNA  release  even  though 
the  macro lide  binding  sites  is  far  below  the  PTC  [26].  It 
has  also  been  reported  that  overexpression  of  very  short 
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Figure  4  Physiological  significance  of  post-initiation  pausing.  (A)  HEK293  cells  expressing  GFP  control,  L4  wild  type,  or 
Aloop  mutant  were  infected  with  lentiviruses  encoding  shRNA  targeting  endogenous  L4  or  scramble  control.  Cell  viability  was 
measured  and  quantified.  Data  are  shown  as  means  ±  SD,  n  =  3.  (B)  A  model  for  post-initiation  pausing  of  ribosomes.  After 
initiation,  the  nascent  chain  (red)  traverses  the  exit  tunnel  (yellow).  At  the  5th  codon  position,  the  paused  ribosome  undergoes 
conformational  changes,  thereby  facilitating  elongation  commitment  by  preventing  nascent  chain  from  off-route  and  abortive 
translation. 


minigenes  triggers  ribosome  drop-off.  However,  when 
the  ORF  length  increased  from  three  codons  to  seven,  the 
drop-off  frequency  decreased  significantly  [27].  These 
findings  can  be  explained  by  the  existence  of  an  addi¬ 
tional  step  between  initiation  and  elongation.  It  will  be 
interesting  to  use  single  molecule  fluorescence  approach 
to  elucidate  ribosome  dynamics  in  real  time  during  the 
first  several  cycles  of  elongation  [28]. 

The  mandatory  pausing  of  ribosomes  at  the  5th  codon 
position  is  advantageous  to  ensure  elongation  commit¬ 
ment.  The  exit  tunnel  is  the  preferred  but  perhaps  not  the 
only  pathway  from  the  PTC.  Some  experimental  results 
suggest  that  the  nascent  chain  could  leave  the  PTC  along 
other  routes  [29].  We  surmise  that  the  post-initiation 
pausing  of  ribosomes  allows  the  very  first  few  amino  ac¬ 
ids  to  And  the  correct  exit  route.  This  unique  interaction 
by  the  first  amino  acids  could  potentially  trigger  local 
conformational  change  of  ribosomes  before  entering  the 
exit  tunnel.  Interestingly,  recent  cryoelectron  microscopy 
of  eukaryotic  ribosomes  revealed  discontinuous  electron 
density  of  nascent  chains  after  ~4  amino  acids  below 
the  PTC  [24],  a  strong  indication  of  ribosome  plasticity. 
In  addition  to  RPL4,  nascent  chains  undergo  extensive 
interaction  with  many  other  ribosomal  proteins,  such  as 
RPL17  and  RPL39.  Biochemical  evidence  suggested  that 
the  tunnel  topology  must  be  flexible  to  accommodate 
the  nascent  chain  [30].  Given  the  unique  feature  of  stalled 
ribosomes  at  the  5th  codon  position,  post-initiation  pausing 
of  ribosomes  represents  a  physiological  signature  of  the 
translation  machinery  to  ensure  productive  translation. 


Materials  and  Methods 

Cells  and  reagents 

HEK293,  HeLa,  and  MEF  cells  were  maintained  in  Dulbecco’s 
Modified  Eagle’s  Medium  (DMEM)  with  10%  fetal  bovine  se¬ 
rum  (FBS).  CHX,  puromycin,  and  erythromycin  were  purchased 
from  Sigma.  LTM  was  provided  by  Ben  Shen  (Scripps  Florida) 
as  described  previously  [19].  Anti-RPS6  antibody  was  purchased 
from  Cell  Signaling  Technology;  anti-myc  from  Santa  Cruz 
Biotechnology;  a  human  polyclonal  autoimmune  antisera  target¬ 
ing  ribosomal  P  is  from  Immunovision;  anti-flag,  anti-RPS20, 
and  anti-RPL5  antibodies  from  Sigma.  Plasmid  transfection  was 
perfonned  using  Lipofectamine  2000  (Invitrogen)  following  the 
manufacturer’s  instruction. 

Plasmids  and  lentiviral  shRNAs 

The  full-length  RPL4  gene  was  amplified  by  PCR  followed  by 
direct  cloning  into  pcDNA3.1-myc-His  vector  (Invitrogen).  To 
construct  the  RPL4(Aloop)  mutant,  RPL4  fragments  correspond¬ 
ing  to  CDS  1-252  nt  and  262-1  281  nt  were  amplified  by  PCR 
followed  by  ligation  using  a  BamHl  restriction  site.  RPL4  shR¬ 
NA  targeting  sequences  (human:  5'-TGAATAAAGACCTGAT- 
TATAC-3')  and  control  sequence  (5'-AACAGTCGCGTTTGC- 
GACTGG-3')  were  designed  based  on  RNAi  consortium  at  Broad 
Institute  (http://www.broad  mit.edu/mai/trc).  All  shRNA  targeting 
sequences  were  cloned  into  DECIPHER™  pRSI9-U6-(sh)-UbiC- 
TagRFP-2A-Puro  (CELLECTA).  Lentiviral  particles  were  pack¬ 
aged  using  Lenti-X  293T  cells  (Clontech).  Vims-containing  su¬ 
pernatants  were  collected  at  48  h  after  transfection  and  target  cells 
were  infected  in  the  presence  of  8  pg/ml  polybrene. 

Immunofluorescence  staining 

Cells  grown  on  glass  coverslips  were  fixed  in  4%  parafonnal- 
dehyde  and  permeabilized  by  0.2%  Triton  X-100.  After  blocking 
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in  2%  BSA  in  PBS,  fixed  cells  were  incubated  with  the  primary 
antibody  at  4  °C  overnight  followed  by  1-h  incubation  at  room 
temperature  with  Alexa  Fluor-labeled  secondary  antibodies.  Cells 
were  then  washed  with  PBS  and  incubated  for  5  min  in  PBS  sup¬ 
plemented  with  Hoechst  to  counter-stain  the  nuclei.  After  a  final 
wash  with  PBS,  coverslips  were  mounted  onto  slides  and  viewed 
using  a  confocal  microscope  (Zeiss  LSM710). 

Immunoprecipitation 

Cells  were  harvested  in  ice-cold  lysis  buffer  (10  mM  Hepes, 
pH  7.4,  100  mM  KCl,  5  mM  MgClj,  100  pg/ml  CHX)  and  lysed 
with  Lysing  Matrix  D  (Fisher).  After  centrifugation  at  4  °C  and 
10  OOOx  g  for  10  min,  supernatant  was  collected,  and  treated  with 
RNase  I  (Ambion)  at  4  °C  for  1  h.  Digestion  was  neutralized  with 
SUPERase  inhibitor  (SUPERase  In,  Ambion)  and  the  samples 
were  incubated  with  anti-myc  beads  (Sigma)  or  anti-flag  beads 
(Sigma)  at  4  °C  overnight.  Immunoprecipitates  were  washed  3 
times  with  lysis  buffer.  The  washed  beads  were  resuspended  in  1  x 
SDS  sample  buffer  (100  mM  Tris,  pH  6.8,  2%  SDS,  15%  glycerol, 
5%  (3-mercaptoethanol,  0.1%  bromophenol  blue),  boiled  for  5  min, 
and  analyzed  by  immunoblotting. 

Immunoblotting 

Proteins  were  resolved  on  SDS-PAGE  and  transferred  to  Immo- 
bilon-P  membranes  (Millipore).  Membranes  were  blocked  for  1  h 
in  TBS  containing  5%  BSA,  followed  by  incubation  with  primary 
antibodies.  After  incubation  with  horseradish  peroxidase-coupled 
secondary  antibodies,  immunoblots  were  developed  using  en¬ 
hanced  chemiluminescence  (ECLPlus,  GE  Healthcare). 

Cell  viability  assay 

HEK293  stable  cell  lines  expressing  GFP  control,  L4  wild  type, 
or  Aloop  mutant  were  infected  with  lentiviruses  encoding  shRNA 
targeting  RPL4  or  scramble  control.  48  h  after  infection,  cells  were 
trypsinized  and  plated  at  equal  density  to  the  96-well  plate.  Cell 
viability  was  measured  on  day  3  using  CCK-8  kit  (Dojindo  Molec¬ 
ular  Technologies)  according  to  the  manufacturer’s  instructions. 

Ribosome  profiling 

Sucrose  solutions  were  prepared  in  polysome  gradient  buffer 
(10  mM  Hepes,  pH  7.4,  100  mM  KCl,  5  mM  MgClj,  100  pg/ml 
CHX,  5  mM  DTT,  and  20  U/ml  SUPERase_In  (Ambion)).  Su¬ 
crose  density  gradients  (15%  -  45%  (wt/vol))  were  freshly  made 
in  SW41  ultracentrifuge  tubes  (Fisher)  using  a  BioComp  Gradient 
Master  (BioComp)  according  to  the  manufacturer’s  instructions. 
Cells  were  plated  to  four  10-cm  dishes  before  ribosome  profiling. 
HEK293  cells  were  first  treated  with  CHX  (100  pM)  or  LTM  (50 
pM)  for  30  min  at  37  °C  to  freeze  the  translating  ribosomes  or  ini¬ 
tiating  ribosomes,  respectively.  After  ice-cold  PBS  solution  wash, 
cells  were  then  harvested  by  ice-cold  polysome  lysis  buffer  (10 
mM  Hepes,  pH  7.4,  100  mM  KCl,  5  mM  MgCl^,  100  pg/ml  CHX, 
5  mM  DTT,  20  U/ml  SUPERaseJn,  and  2%  (vol/vol)  Triton 
X-100).  After  centrifugation  at  4  °C  and  10  OOOx  g  for  10  min,  ap¬ 
proximately  650  pi  supernatant  was  loaded  onto  sucrose  gradients, 
followed  by  centrifugation  for  100  min  at  38  000  rpm,  4  °C,  in 
an  SW41  rotor.  Separated  samples  were  fractionated  at  0.375  ml/ 
min  by  using  a  fractionation  system  (Isco)  that  continually  moni¬ 
tored  OD254  values.  Fractions  were  collected  into  tubes  at  1-min 
intervals.  To  convert  the  polysome  into  monosome,  E.  coli  RNase 


I  (Ambion)  was  added  into  the  pooled  polysome  samples  (750  U 
per  100  A260  units)  and  incubated  at  4  °C  for  1  h.  SUPERase_In 
(50  U  per  100  U  RNase  I)  was  then  added  to  stop  digestion.  Total 
RNA  extraction  was  perfonned  using  TRIzol  reagent. 

cDNA  library  construction  of  RPF 

Purified  RNA  samples  were  dephosphorylated  in  a  15  pi  reac¬ 
tion  containing  lx  T4  polynucleotide  kinase  buffer,  10  U  SUPER- 
ase_In,  and  20  U  T4  polynucleotide  kinase  (NEB).  Dephosphor¬ 
ylation  was  carried  out  for  1  h  at  37  °C,  and  the  enzyme  was  then 
heat  inactivated  for  20  min  at  65  °C.  Dephosphorylated  samples 
were  then  mixed  with  2  Novex  TBE-Urea  sample  buffer  (Invitro- 
gen)  and  loaded  on  a  Novex  denaturing  15%  polyacrylamide  TBE- 
urea  gel  (Invitrogen).  The  gel  was  stained  with  SYBR  Gold  (Invit- 
rogen)  to  visualize  the  RNA  fragments.  Gel  bands  containing  RNA 
species  corresponding  to  28  nt  were  excised  and  physically  dis¬ 
rupted  by  centrifugation  through  the  holes  of  the  tube.  RNA  frag¬ 
ments  were  dissolved  by  soaking  overnight  in  gel  elution  buffer 
(300  mM  NaOAc,  pH  5.5,  1  mM  EDTA,  0.1  U/ml  SUPERaseJn). 
The  gel  debris  was  removed  using  a  Spin-X  column  (Corning)  and 
RNA  was  purified  by  using  ethanol  precipitation.  Purified  RNA 
fragments  were  resuspended  in  10  mM  Tris  (pH  8)  and  denatured 
briefiy  at  65  °C  for  30  s.  Poly-(A)  tailing  reaction  was  performed 
in  a  8  pi  buffer  with  lx  poly-(A)  polymerase,  1  mM  ATP,  0.75  U/ 
pi  SUPERaseJn,  and  3  U  E.  coli  poly-(A)  polymerase  (NEB). 
Tailing  was  carried  out  for  45  min  at  37  °C.  For  reverse  transcrip¬ 
tion,  the  following  oligos  containing  barcodes  were  synthesized: 

MCA02,  5'-pCAGATCGTCGGACTGTAGAACTCT/idSp/CAAG- 
CAGAAGACGGCATACGATTTTTTTTTTTTTTTTTTTTVN-3'; 
LGT03,  5'-pGTGATCGTCGGACTGTAGAACTCT/idSp/CAAG- 
CAGAAGACGGCATACGATTTTTTTTTTTTTTTTTTTTVN-3'; 
YAG04,  5  '-pTCGATCGTCGGACTGTAGAACTCT/idSp/C AAG- 
CAGAAGACGGCATACGATTTTTTTTTTTTTTTTTTTTVN-3'; 
HTC05,  5'-pAGGATCGTCGGACTGTAGAACTCT/idSp/CAAG- 
CAGAAGACGGCATACGATTTTTTTTTTTTTTTTTTTTVN-3'. 

In  brief,  the  tailed  RNA  product  was  mixed  with  0.5  mM  dNTP 
and  2.5  mM  synthesized  primer  and  incubated  at  65  °C  for  5  min, 
followed  by  incubation  on  ice  for  5  min.  The  reaction  mix  was 
then  added  with  20  mM  Tris  (pH  8.4),  50  mM  KCl,  5  mM  MgClj, 
10  mM  DTT,  40  U  RNaseOUT,  and  200  U  Superscript  III  (Invit¬ 
rogen).  Reverse  transcription  reaction  was  perfonned  according  to 
the  manufacturer’s  instructions.  RNA  was  eliminated  from  cDNA 
by  adding  1 .8  pi  of  1  M  NaOH  and  incubating  at  98  °C  for  20  min. 
The  reaction  was  then  neutralized  with  1 .8  pi  of  1  M  HCl.  Reverse 
transcription  products  were  separated  on  a  10%  polyacrylamide 
TBE-urea  gel  as  described  earlier.  The  extended  first-strand  prod¬ 
uct  band  was  expected  to  be  approximately  100  nt,  and  the  corre¬ 
sponding  region  was  excised.  The  cDNA  was  recovered  by  using 
DNAgel  elution  buffer  (300  mM  NaCl,  1  mM  EDTA).  First-strand 
cDNA  was  circularized  in  20  pi  of  reaction  containing  1  x  CircLi- 
gase  buffer,  2.5  mM  MnC12,  IM  Betaine,  and  100  U  CircLigase 

II  (Epicentre).  Circularization  was  performed  at  60  °C  for  1  h,  and 
the  reaction  was  heat  inactivated  at  80  °C  for  10  min.  Circular  sin¬ 
gle-strand  DNA  was  relinearized  with  20  mM  Tris-acetate,  50  mM 
potassium  acetate,  10  mM  magnesium  acetate,  1  mM  DTT,  and 
7.5  U  APE  1  (NEB).  The  reaction  was  carried  out  at  37  °C  for  1  h. 
The  linearized  single-strand  DNA  was  separated  on  a  Novex  10% 
polyacrylamide  TBE-urea  gel  (Invitrogen)  as  described  earlier. 
The  expected  100-nt  product  bands  were  excised  and  recovered  as 
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described  earlier. 

Deep  sequencing 

Single-stranded  template  was  amplified  by  PCR  using  the  Phu- 
sion  High-Fidelity  enzyme  (NEB)  according  to  the  manufacturer’s 
instructions.  The  oligonucleotide  primers  qNTI200  (5'-CAAG- 
CAGAAGACGGCATA-3')  and  qNTI201  (5'-AATGATACGGC- 
GACCACCGACAGGTTCAGAGTTCTACAGTCCGACG-3') 
were  used  to  create  DNA  suitable  for  sequencing,  i.e.,  DNA  with 
Illumina  cluster  generation  sequences  on  each  end  and  a  sequenc¬ 
ing  primer  binding  site.  The  PCR  contains  lx  HF  buffer,  0.2  mM 
dNTP,  0.5  pM  oligonucleotide  primers,  and  0.5  U  Phusion  poly¬ 
merase.  PCR  was  carried  out  with  an  initial  30-s  denaturation  at 
98  °C,  followed  by  12  cycles  of  10-s  denaturation  at  98  °C,  20-s 
annealing  at  60  °C,  and  10  s  extension  at  72  °C.  PCR  products 
were  separated  on  a  nondenaturing  8%  polyacrylamide  TBE  gel 
as  described  earlier.  Expected  DNA  at  120  bp  was  excised  and 
recovered  as  described  earlier.  After  quantification  by  Agilent 
BioAnalyzer  DNA  1000  assay,  equal  amount  of  barcoded  samples 
were  pooled  into  one  sample.  Approximately  3-5  pM  mixed  DNA 
samples  were  used  for  cluster  generation  followed  by  sequencing 
using  sequencing  primer  5'-CGACAGGTTCAGAGTTCTA- 
CAGTCCGACGATC-3'  (Illumina  Genome  Analyzer  2  or  HiSeq). 

Data  analysis 

Barcoded  reads  were  separated  into  individual  data  file  for  each 
sample  according  to  the  5'  end  2-nt  barcode  sequences.  Seven 
nucleotides  were  trimmed  from  the  3'  end  of  each  50-nt-long  Illu¬ 
mina  sequence  read  to  remove  adaptor  sequences,  and  the  poly  (A) 
tails  were  removed  from  the  3'  end,  allowing  one  mismatch.  Reads 
between  25  and  35  nt  in  length  were  mapped  using  Bowtie  [31]  to 
the  sense  strand  of  the  longest  Refseq  transcripts  for  each  human 
gene.  One  mismatch  was  allowed  in  all  mappings;  in  cases  of 
multiple  mapping,  mismatched  positions  were  not  used  if  a  perfect 
match  existed.  Reads  mapped  more  than  100  times  were  discarded 
to  remove  poly-A-derived  reads.  Finally,  reads  were  counted  at  ev¬ 
ery  position  of  individual  transcripts  by  using  the  13th  nucleotide 
of  the  read  for  the  P-site  position. 

To  generate  aggregation  plots,  the  number  of  RPF  reads  aligned 
to  each  position  of  individual  transcripts  was  first  normalized  by 
the  read  density  within  50-150  codon  window,  which  was  consid¬ 
ered  as  the  representative  translational  level  of  the  transcript.  Only 
transcripts  with  CDS  longer  than  150  codons  and  with  minimum 
read  density  of  25  rpkM  were  used  for  aggregation  plots.  The  read 
counts  were  then  averaged  across  all  transcripts  for  each  position 
relative  to  the  annotated  start  codon  based  on  Refseq. 

The  framing  of  reads  mapped  to  the  1st  codon,  5th  codon,  or 
entire  CDS  was  calculated  relative  to  the  annotated  ORF  based 
on  Refseq.  The  length  distribution  was  also  calculated  for  reads 
mapped  to  the  1st  codon,  the  5th  codon,  or  the  entire  CDS.  Read 
density  over  the  CDS  (6-stop  codon)  relative  to  the  initiation  site 
(1-5  codon)  was  computed  for  each  transcript.  Only  transcripts 
with  minimum  read  density  of  10  rpkM  over  the  CDS  and  mini¬ 
mum  5  reads  at  the  initiation  site  were  used  to  calculate  the  ratios. 

For  codon  composition  analysis,  the  annotated  CDS  from  hu¬ 
man  gene  was  downloaded  from  UCSC  table  browser.  The  codon 
composition  of  the  1st  codon,  5th  codon,  and  whole  CDS  were 
summarized  based  on  the  longest  transcripts  of  human  genes. 

For  GO  terai  enrichment  analysis,  two  groups  of  genes  were 
collected  based  on  the  average  read  densities  in  the  regions  Hank¬ 


ing  the  start  codon  and  the  5th  codon  by  the  following  criteria: 

Group  I:  Density >  1.2  x  Density 

Group  II:  Density >  1.2  x  Density 

The  gene  ontology  analyses  were  conducted  on  the  above  two 
groups  of  genes  using  DAVID  [32]. 
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Cells  have  evolved  exquisite  mechanisms  to  fine-tune  the 
rate  of  protein  synthesis  in  response  to  stress.  Systemic 
mapping  of  start-codon  positions  and  precise  measurement 
of  the  corresponding  initiation  rate  would  transform  our 
understanding  of  translational  control.  Here  we  present 
quantitative  translation  initiation  sequencing  (QTI-seq), 
with  which  the  initiating  ribosomes  can  be  profiled  in  real  time 
at  single-nucleotide  resolution.  Resultant  initiation  maps  not 
only  delineated  variations  of  start-codon  selection  but  also 
highlighted  a  dynamic  range  of  initiation  rates  in  response 
to  nutrient  starvation.  The  integrated  data  set  provided 
unique  insights  into  principles  of  alternative  translation 
and  mechanisms  controlling  different  aspects  of  translation 
initiation.  With  RiboTag  mice,  QTI-seq  permitted  tissue- 
specific  profiling  of  initiating  ribosomes  in  vivo.  Liver 
cell-specific  ribosome  profiling  uncovered  a  robust 
translational  reprogramming  of  the  proteasome  system 
in  fasted  mice.  Our  findings  illuminated  the  prevalence 
and  dynamic  nature  of  translational  regulation  pivotal  to 
physiological  adaptation  in  vivo. 

Translational  regulation  permits  cells  to  respond  swiftly  to 
stress  conditions  via  immediate  and  selective  changes  in  protein 
expression  levels^”^.  Much  of  this  translational  control  occurs 
at  the  initiation  stage^“^.  Although  the  importance  of  transla¬ 
tion  initiation  has  long  been  appreciated,  direct  measurement 
of  the  initiation  rate  on  individual  mRNAs  has  been  difficult. 
Recent  development  of  ribosome-profiling  technology  (Ribo- 
seq),  based  on  deep  sequencing  of  ribosome-protected  mRNA 
fragments  (RPFs),  enables  monitoring  of  ribosome  dynamics 
with  unprecedented  resolution  at  the  genome-wide  scale®’®.  The 
obtained  snapshot  of  ribosome  occupancy  on  the  coding  region 
has  often  been  used  to  estimate  relative  changes  of  translation 
efficiency  under  different  growth  conditions.  A  caveat,  however, 
exists  because  the  average  ribosome  density  on  mRNAs  is  nega¬ 
tively  influenced  by  the  elongation  speed^®’^*.  Dynamic  changes 
of  initiation  rates  are  thus  masked  by  the  varied  elongation  speed 
under  different  growth  conditions.  In  addition,  the  average  den¬ 
sity  of  ribosome  occupancy  is  not  suitable  to  evaluate  alterna¬ 
tive  translation  events  occurring  on  the  same  mRNA^^.  These 
constraints  point  to  a  need  for  a  method  capable  of  mapping 


start-codon  selection  and  quantifying  the  rate  of  SOS  ribosome 
assembly  at  individual  translation  initiation  sites  (TISs). 

Much  of  our  understanding  of  the  basic  mechanisms  of 
translation  initiation  is  limited  to  cells  in  culture.  Investigating 
translational  control  at  the  organismal  level  remains  a  pressing 
and  formidable  challenge.  The  cellular  heterogeneity  of  tissues 
and  organs  confounds  our  efforts  to  achieve  cell  type-specific 
genomic  and  proteomic  interrogation.  To  monitor  protein  syn¬ 
thesis  at  the  organismal  level,  one  approach  relies  on  tagged 
ribosomal  proteins,  such  as  hemagglutinin  (HA) -tagged  Rpl22 
(ref  13)  or  EGFP-fused  RpllOa  (ref  14).  Following  genetic  target¬ 
ing  to  specific  cell  populations,  affinity  purification  allows  isola¬ 
tion  of  polysome-bound  mRNAs  from  specific  cell  types  suitable 
for  profiling  assays.  However,  current  methods  do  not  allow 
tissue-specific  profiling  of  initiating  ribosomes,  which  bear  more 
reliable  information  in  evaluating  translational  control  in  vivo. 

Here  we  present  an  approach  called  QTI-seq  that  allows 
comprehensive  monitoring  of  translation  initiation  in  cells  and 
solid  tissues.  Together  with  the  RNA-seq  and  Ribo-seq  data  sets 
acquired  in  parallel,  the  integrated  ribosome-profiling  data  afford 
a  holistic  view  of  the  translational  landscape  under  different 
growth  conditions. 

RESULTS 

QTI-seq  captures  initiating  ribosomes  in  an  unbiased  manner 

We  previously  developed  global  translation  initiation  sequencing 
(GTI-seq)^®,  which  capitalizes  on  the  unique  feature  of  lactimid- 
omycin  (LTM)  that  preferentially  acts  on  the  first  SOS  ribosome 
when  its  E  site  is  empty^®.  However,  LTM-mediated  enrichment 
of  initiating  ribosomes  requires  an  incubation  of  the  cell  at  37  °C 
to  allow  for  the  release  of  elongating  ribosomes.  This  procedure, 
albeit  essential  for  high-resolution  mapping  of  TISs,  may  gener¬ 
ate  a  number  of  artifacts.  First,  free  ribosomes  could  start  new 
rounds  of  initiation  during  the  incubation  period,  resulting 
in  an  amplification  of  ribosome  occupancy  at  the  start  codon 
(Supplementary  Fig.  la).  Second,  the  amplified  ribosome  density 
is  biased  toward  the  upstream  initiators,  causing  an  unwanted 
5'-end  inflation  of  RPFs.  When  GTI-seq  was  applied  to  human 
HEK293  cells  with  amino  acid  deprivation,  the  LTM  peaks  at  the 
annotated  start  codons  (aTISs)  showed  few  changes  despite  the 
pervasive  reduction  of  ribosome  occupancy  on  many  transcripts 
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Figure  1  |  QTI-seq  captures  real-time  translation  initiation  events  in  a  qualitative  and  quantitative  manner,  (a)  Schematic  of  Ribo-seq  (left)  and 
QTI-seq  (right)  procedures.  CHX,  cycloheximide;  LTM,  lactimidomydn;  PMY,  puromycin.  (b)  Meta-gene  analysis  of  CHX-associated  ribosome  density 
and  LTM-assodated  ribosome  density  in  HEK293  cells  captured  by  QTI-seq.  Normalized  ribosome-protected  mRNA  fragment  (RPF)  reads  are  averaged 
across  the  entire  transcriptome  and  aligned  at  the  annotated  start  codon  (aTIS).  Different  reading  frames  are  separated  and  color  coded,  (c)  Relative 
ratios  of  different  types  of  TISs  identified  by  GTI-seq  or  QTI-seq  in  HEK293  cells.  uTIS,  upstream  TIS;  dTIS,  downstream  TIS.  (d)  Correlation  between 
LTM-assodated  aTIS  density,  normalized  by  upper  quartile  (UQ),  and  CHX-assodated  coding  region  (CDS)  ribosome  occupancy,  normalized  by  reads  per 
kilobase  per  million  reads  (RPKM).  Genes  with  a  single  aTIS  or  multiple  TISs  are  shown  as  blue  and  red  dots,  respectively,  (e)  Examples  of  single-TIS 
{ABCD3)  and  multiple-TIS  {6EMIN6)  genes  revealed  by  GTI-seq  (top)  and  QTI-seq  (bottom).  They-axis  scales  are  identicaL  The  corresponding  gene 
structure  is  shown  below  the  x  axis. 


(Supplementary  Fig.  lb).  As  a  result,  virtually  no  correlation 
was  found  between  the  differences  of  LTM-captured  aTIS  den¬ 
sity  and  the  ribosome  occupancy  on  the  coding  region  (CDS) 
(Supplementary  Fig.  Ic).  The  lack  of  quantitative  features  in 
GTI-seq  largely  hampers  comparative  interrogation  of  transla¬ 
tional  regulation. 

To  circumvent  this  limitation,  we  set  out  to  devise  QTI-seq 
to  preserve  the  small  population  of  initiating  ribosomes  with 
minimal  perturbation.  Following  extensive  optimization,  we 
implemented  three  steps  to  achieve  this  goal:  (i)  rapidly  break 
down  cells  using  lysing  matrix  D,  which  has  a  minimal  effect  on 
ribosome  stability;  (ii)  specifically  freeze  the  initiating  ribosomes 
by  treating  cell  lysates  with  LTM;  and  (iii)  effectively  deplete  the 
elongating  ribosomes  by  introducing  another  translation  inhibi¬ 
tor,  puromycin  (PMY)  (Fig.  la).  PMY  does  not  affect  translation 
initiation  and  has  been  used  previously  to  enrich  the  initiating 
ribosomes*^.  Acting  as  a  tRNA  analog,  PMY  releases  the  nas¬ 
cent  chain  and  dissociates  the  ribosome  into  separate  subunits’®. 
In  the  presence  of  elongation  inhibitors  such  as  cycloheximide 
(CHX),  however,  addition  of  PMY  catalyzes  puromycylation  of 
the  nascent  chain  without  triggering  ribosome  dissociation’®. 
The  same  feature  holds  true  for  cell-free  samples  as  the  pres¬ 
ence  of  CHX  completely  prevented  PMY- induced  polysome  dis¬ 
assembly  (Supplementary  Fig.  2).  Because  LTM  uses  a  similar 
mechanism  to  that  of  CHX  but  selectively  acts  on  the  first  SOS 
ribosome’®,  we  reasoned  that  sequential  treatment  with  LTM  and 
PMY  would  dissociate  the  elongating  ribosomes  but  leave  the 
initiating  ribosomes  insensitive  to  PMY  (Fig.  la).  Indeed,  PMY 


addition  in  the  presence  of  LTM  resulted  in  polysome  disassembly 
with  a  corresponding  increase  of  the  monosome  (Supplementary 
Fig.  2).  To  examine  whether  the  increased  monosome  contained 
the  preserved  initiating  ribosomes,  we  purified  RPFs  followed 
by  deep  sequencing.  Metagene  analysis  revealed  that  the  LTM- 
associated  RPF  reads  were  highly  enriched  at  the  aTIS  codon 
(Fig.  lb).  Additionally,  the  aTIS  ribosome  density  captured  by 
LTM  was  highly  reproducible  (Supplementary  Fig.  3). 

Because  QTl-seq  captured  a  small  population  of  initiating 
ribosomes  without  5' -end  RPF  inflation,  QTI-seq  identified  fewer 
total  TISs  than  GTI-seq  in  HEK293  cells  (5,099  vs.  13,915),  espe¬ 
cially  the  upstream  TIS  codons  (uTISs)  (Fig.  Ic).  Consistent  with 
previous  reports’®'^®,  codon  composition  analysis  revealed  that 
more  than  half  the  TISs  used  a  non- AUG  sequence  as  the  transla¬ 
tion  initiator  (Supplementary  Fig.  4).  The  associated  LTM  peak 
height  at  the  aTISs  showed  a  strong  correlation  with  the  average 
ribosome  occupancy  along  the  corresponding  CDS  (r  =  0.728; 
Fig.  Id).  In  addition  to  this  quantitative  feature,  QTI-seq  retained 
the  high  precision  in  mapping  TIS  positions  at  a  single-nucleotide 
resolution.  For  instance,  a  prominent  LTM  peak  was  located 
exactly  at  the  annotated  start  codon  of  ABCD3;  an  out-of- frame 
uTIS  was  uncovered  in  GEMIN6  (Fig.  le).  QTI-seq  thus  offers  a 
promising  approach  to  exploring  real-time  translation  initiation 
in  a  qualitative  and  quantitative  manner. 

Quantitative  TIS  profile  in  response  to  starvation 

We  next  applied  QTI-seq  to  HEK293  cells  with  amino  acid  dep¬ 
rivation  (Supplementary  Fig.  5).  Total  cellular  RNA  was  also 
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Figure  2  \  QTI-seq  reveals  pervasive  translational  regulation  in 
response  to  starvation.  (a,b)  Fold  changes  in  LTM-associated  aTIS 
density  and  CHX-associated  CDS  ribosome  occupancy  in  HEK293  cells  (a) 
and  MEFs  (b)  before  and  after  amino  add  starvation.  Genes  with 
more  than  1.5-fold  changes  are  within  "Increase"  or  "Decrease"  boxes. 

(c)  Fleat  map  of  false  discovery  rate  (FDR)  of  enriched  Gene  Ontology 
(GO)  "Biological  process"  terms  for  gene  groups  with  translational 
downregulation  (Dec)  or  upregulation  (Inc)  in  response  to  amino  add 
starvation,  (d)  Left,  example  of  a  single-TIS  gene  (RPS28)  showing 
a  decreased  aTIS  after  amino  acid  starvation.  Right,  relative  RNA 
abundance,  CDS  ribosome  occupancy  and  aTIS  ribosome  density  of 

RPS28.  (e)  Example  of  a  multiple-TIS  gene  {6ADD456)  showing  an  increased  aTIS  fraction  after  amino 
translational  control  of  GADD45G  by  a  firefly  luciferase  (Flue)  reporter  bearing  the  5'  UTR  of  GADD45G 
Data  are  mean  ±  s.e.m.;  n  «  3  biological  replicates;  *P  <  0.05  by  Student's  t-test. 
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acid  starvation,  (f)  Experimental  validation  of 
with  (WT)  or  without  the  uTIS  codon  CTG  (Del). 


collected  in  parallel  for  RNA-seq  to  quantify  mRNA  abundance. 
In  response  to  starvation,  the  changes  in  the  initiation  rates 
showed  positive  correlation  with  the  differences  of  ribosome 
occupancy  on  the  corresponding  CDS  (r  =  0.375;  Fig.  2a).  The 
same  observation  held  for  a  mouse  embryonic  fibroblast  (MEF) 
cell  line  subjected  to  starvation  (r  =  0.419;  Fig.  2b).  The  imperfect 
correlation  was  partially  due  to  reduced  elongation  speed  under 
nutrient  starvation  that  potentially  increases  the  CDS  ribosome 
occupancy  (Supplementary  Fig.  6).  From  the  comprehensive  data 
sets  acquired  from  QTI-seq  and  Ribo-seq,  we  identified  a  large 
number  of  transcripts  that  underwent  twofold  changes  upon  star¬ 
vation  (1,073  in  HEK293  cells  and  820  in  MEFs;  Supplementary 
Tables  1  and  2).  Among  the  genes  showing  repressed  transla¬ 
tion,  many  of  them  are  involved  in  protein  biosynthesis  and 
metabolism  (Fig.  2c).  As  a  typical  example,  the  gene  encoding 
ribosomal  protein  RPS28  showed  a  nearly  fivefold  decrease 
in  ribosome  occupancy  on  the  CDS  in  response  to  starvation 
(Fig.  2d).  Remarkably,  QTI-seq  displayed  a  decrease  greater 
than  14-fold  in  the  ribosome  density  at  the  start  codon  of  RPS28. 
Without  the  negative  influence  imposed  by  elongation  speed  vari¬ 
ation,  QTI-seq  sets  itself  apart  from  previous  methods  by  report¬ 
ing  true  measures  of  initiation  rates. 

QTI-seq  also  captured  a  considerable  number  of  genes  whose 
transcripts  experienced  increased  initiation  in  response  to  star¬ 
vation  (638  in  HEK293  cells  and  515  in  MEFs).  Interestingly, 
many  of  these  genes  are  involved  in  protein  folding  and  nuclear 
transport  (Fig.  2c).  Using  a  firefly  luciferase  (Flue)  reporter  con¬ 
taining  the  5'  UTR,  we  validated  the  translational  upregulation 
of  nucleoporin-encoding  gene  NUP88  (Supplementary  Fig.  7). 
It  is  noteworthy  that  many  starvation-responsive  genes  contain 
multiple  TISs  (1,286  in  HEK293  cells  and  1,343  in  MEFs),  which 
suggests  a  regulatory  role  for  alternative  TISs  in  translational 
controFh  To  demonstrate  the  influence  of  alternative  translation 


on  the  aTIS  initiation,  we  selected  genes  with  multiple  TISs  and 
computed  the  ribosome  density  at  the  aTIS  codon  relative  to  the 
total  TIS  density  on  the  same  transcript  (Supplementary  Fig.  8). 
This  analysis  uncovered  many  genes  whose  translational  regula¬ 
tion  is  indiscernible  by  simple  comparison  of  ribosome  density 
changes  at  either  CDSs  or  aTISs.  A  total  of  428  genes  in  HEK293 
cells  and  212  genes  in  MEFs  demonstrated  an  altered  aTIS-to-TIS 
ratio  upon  amino  acid  deprivation  (false  discovery  rate  <0.05). 
This  strategy  uncovered  several  stress-responsive  genes  whose 
transcripts  contain  previously  uncharacterized  TISs.  For  instance, 
GADD45G  bears  a  CUG  start  codon  in  the  5'  UTR  (Fig.  2e).  Flue 
reporter  assays  confirmed  the  critical  role  of  5'  UTR  in  the  starva¬ 
tion-induced  upregulation  of  GADD45G  (Fig.  2f).  In  particular, 
deleting  the  CUG  codon  was  sufficient  to  prevent  the  starvation 
responsiveness. 

Programmatic  TIS  regulation  in  response  to  starvation 

Many  upstream  open  reading  frames  are  believed  to  exert  nega¬ 
tive  effects  on  the  main  open  reading  frame  translation,  pre¬ 
sumably  by  capturing  the  scanning  ribosome It  is  thus  not 
surprising  to  find  that  a  large  number  of  multiple- TIS-containing 
genes  showed  increased  aTIS  initiation  when  uTIS  initiation  was 
repressed  under  starvation.  However,  a  handful  of  transcripts 
exhibited  decreased  aTIS  fraction  despite  the  presence  of  alter¬ 
native  TISs  (Supplementary  Fig.  8).  To  identify  possible  factors 
governing  differential  regulation  of  alternative  TISs,  we  surveyed 
for  consensus  sequence  motifs  in  gene  groups  that  respond  differ¬ 
ently  to  starvation.  Among  transcripts  with  increased  aTIS  initia¬ 
tion  upon  starvation,  the  Kozak  consensus  motif  was  prominent 
(Supplementary  Fig.  9a).  For  transcripts  with  increased  uTISs, 
an  evident  purine-rich  sequence  context  emerged  above  the  back¬ 
ground  (Supplementary  Fig.  9b).  This  finding  is  reminiscent  of 
the  polypurine  ( A)-rich  sequences  found  in  many  evolutionarily 
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Figure  3  |  Distinct  role  of  eIF2a  phosphorylation  in  translational  response  to  starvation,  (a)  Left,  schematic  of  mTORCl  signaling  and  eIF2oc 
phosphorylation  response.  Right,  immunoblotting  results  of  eIF2a(S51D)  mutant  cells  with  either  amino  acid  starvation  or  doxycycline  (DOX)-induced 
eIF2a(S51D)  expression,  (b)  Fold  changes  in  LTM-associated  aTIS  density  in  eIF2a(S51D)  mutant  cells  between  amino  add  starvation  and  DOX-induced 
eIF2(x(S51D)  expression.  Terminal  oligopyrimidine  (TOP)  mRNAs  are  shown  as  red  dots,  (c)  Heat  map  of  fold  change  of  enriched  GO  terms 
("Biological  process")  for  gene  groups  with  translational  downregulation  (blue)  or  upregulation  (red)  in  response  to  amino  add  starvation  or 
DOX-induced  eIF2a(S51D)  expression,  (d)  Left,  example  of  a  gene  (RPSll)  showing  translational  repression  in  response  to  amino  acid  starvation 
but  not  DOX-induced  eIF2a(S51D)  expression.  Right,  validation  of  RPSll  translational  control  by  a  Flue  reporter  bearing  the  5'  UTR  of  RPSll. 

(e)  Left,  example  of  a  gene  {NUP43)  showing  translational  upregulation  in  response  to  either  amino  acid  starvation  or  DOX-induced  eIF2a(S51D) 
expression.  Right,  validation  of  NUP43  translational  control  by  a  Flue  reporter  bearing  the  5'  UTR  of  NUP43.  Data  in  d,e  are  mean  ±  s.e.m.; 
n  -  3  biological  replicates;  *P  <  0.01  by  Student's  t-test. 


conserved  internal  ribosome  entry  sltes^^'^^.  Alternatively,  the 
purine-rich  sequence  may  act  in  a  manner  opposite  to  the  5'  ter¬ 
minal  oligopyrimidine  (TOP)  tract,  whose  translation  is  highly 
sensitive  to  the  mTORCl  signaling  pathway^^.^s. 

TIS  regulation  pathways  in  response  to  starvation 
Amino  acid  deprivation  suppresses  global  protein  synthesis 
by  inhibiting  the  mTORCl  signaling  pathway  and  activating 
GCN2  kinases^^.  Whereas  the  former  regulates  cap  recognition, 
the  latter  controls  the  ternary  complex  formation  (Fig.  3a).  It 
is  unclear  whether  these  distinct  signaling  pathways  have  over¬ 
lapping  or  unique  contributions  to  the  translational  regulation. 
We  took  advantage  of  a  HEK293  cell  line  harboring  an  inducible 
phosphor- mimetic  allele  ofeIF2a  in  which  serine  51  was  mutated 
to  an  aspartic  acid  (SSID)^®.  Expression  of  this  mutant  upon  doxy¬ 
cycline  (DOX)  addition  led  to  translational  attenuation  as  well 
as  an  evident  induction  of  ATF4  (Fig.  3a  and  Supplementary 
Fig.  10a).  Unlike  amino  acid  deprivation,  however,  eIF2a(S51D) 
expression  had  little  effect  on  mTORCl  signaling  (Fig.  3a). 
Therefore,  the  direct  contribution  of  eIF2a  phosphorylation  to 
translational  reprogramming  can  be  dissected  without  pleiotropic 
effects  caused  by  nutrient  deprivation. 

We  applied  QTI-seq  to  HEK293/eIF2a(S5lD)  cells  with  or 
without  DOX  treatment,  together  with  RNA-seq  and  Ribo-seq. 
Similarly  to  amino  acid  starvation,  eIF2a(S5lD)  expression 


resulted  in  evident  differences  for  the  LTM-associated  aTIS 
ribosome  density  captured  by  QTI-seq  (Supplementary 
Fig.  10b  and  Supplementary  Table  3).  A  comparative  analysis  of 
the  aTIS  profiles  between  these  two  conditions  revealed  a  partially 
overlapping  gene  pattern  (Fig.  3b).  Interestingly,  the  majority  of 
overlapping  genes  were  translationally  upregulated  with  increased 
aTIS  ribosome  density.  For  those  mRNAs  undergoing  transla¬ 
tional  repression  under  starvation,  very  few  were  affected  in  the 
presence  of  eIF2a(S5lD),  with  the  most  notable  gene  group  as 
TOP  mRNAs  (Fig.  3b).  Clustering  of  the  transcripts  with  changed 
aTIS  ribosome  density  further  confirmed  the  dichotomy  of  regu¬ 
lation  between  eIF2a(S5lD)  expression  and  starvation  (Fig.  3c). 
In  general,  the  transcripts  commonly  affected  under  both  condi¬ 
tions  were  translationally  upregulated  and  span  diverse  functions, 
whereas  the  starvation-sensitive  mRNA  subsets  were  translation- 
ally  repressed  and  enriched  in  functions  highly  relevant  to  ribo¬ 
some  biogenesis.  These  results  indicate  that,  in  response  to  amino 
acid  starvation,  reduced  mTORCl  signaling  is  more  responsible 
for  the  suppression  of  TOP  mRNA  translation.  In  contrast,  eIF2a 
phosphorylation  contributes  mainly  to  the  selective  upregulation 
of  a  subset  of  transcripts. 

To  confirm  the  differential  regulation  of  translation  initiation 
controlled  by  distinct  upstream  signaling  pathways,  we  chose  two 
genes  for  experimental  validation.  The  gene  encoding  ribosomal 
protein  RPSll  contains  a  typical  TOP  sequence  on  its  5'  UTR. 
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Figure  4  |  Liver-spedfic  QTI-seq  reveals  translational  reprogramming  in 
response  to  fasting,  (a)  Metagene  analysis  of  LTM-assodated  ribosome 
density  (top)  or  CHX-assodated  ribosome  density  (bottom)  in  MEFs 
and  liver  cells.  Normalized  RPF  reads  are  averaged  across  the  entire 
transcriptome  and  aligned  at  the  annotated  start  codons  and  stop 
codons,  (b)  Fold  changes  in  LTM-associated  aTIS  density  and  CHX- 
associated  CDS  ribosome  occupancy  in  mouse  liver  cells  with  and  without 
fasting,  (c)  Heat  map  of  fold  changes  for  gene  groups  with  translational 
downregulation  (blue)  or  upregulation  (red)  in  liver  cells  of  fasted  mice 
or  starved  MEFs.  (d)  Reporter  assay  using  an  in  vitro  translation  system 
reprogrammed  from  mouse  liver  lysates  with  or  without  fasting. 

The  relative  translation  efficiency  of  a  synthesized  Flue  mRNA  containing 
5'  UTRs  of  Psma3  or  Psmb4  is  shown  along  with  a  no-5'  UTR  condition, 
(e)  Chymotrypsin  activity  of  liver  homogenates  from  mice  8-12  weeks 
old  that  were  either  fed  or  food  deprived  overnight,  measured  by 
Proteasome-Glo.  Data  in  d,e  are  mean  ±  s.e.m.;  n  =  3  biological 
replicates;  *P  <  0.01  by  Student's  t-test. 


Amino  acid  starvation  led  to  an  approximately  threefold  decrease 
of  the  aTIS  ribosome  density  (Fig.  3d).  However,  the  translation 
of  RPSll  exhibited  little  change  in  response  to  DOX-induced 
eIF2a(S5lD)  expression.  The  Flue  reporter  bearing  the  RPSll  5' 
UTR  behaved  exactly  like  the  endogenous  RPSl  i  in  transfected 
cells  (Fig.  3d).  Amino  acid  starvation,  but  not  eIF2a(S5lD) 
expression,  caused  a  significant  reduction  of  Flue  expression. 
Unlike  RPSll,  the  nucleoporin  gene  NUP43  exhibited  a  transla¬ 
tional  upregulation  after  starvation  or  eIF2a(S5  ID)  expression  as 
evidenced  by  a  more-than-threefold  increase  for  its  aTIS  ribosome 
density  (Fig.  3e).  In  agreement  with  the  QTI-seq  data,  the  Flue 
reporter  fused  with  the  5'  UTR  of  NUP43  showed  a  consistent 
increase  of  translation  in  response  to  both  conditions  (Fig.  3e). 
These  results  strongly  favor  the  hypothesis  that  distinct  upstream 
signaling  pathways  act  on  different  aspects  of  translation  initia¬ 
tion,  resulting  in  a  coordinated  translational  reprogramming  to 
achieve  cellular  adaptation. 

Tissue-specific  QTI-seq  in  liver  cells 

Having  elucidated  translational  regulation  from  cells  in  culture, 
we  reasoned  that  the  principle  of  QTI-seq  is  applicable  to  solid  tis¬ 
sues  because  it  directly  captures  initiating  ribosomes  from  lysates. 
To  achieve  cell  type-specific  ribosome  profiling,  we  made  use 
of  the  available  RiboTag  mouse,  which  carries  a  ‘floxed’  Rpl22^ 
allele^^.  As  a  pilot  study,  we  crossed  the  RiboTag  mouse  with  a 
mouse  expressing  liver  cell-specific  Cre  recombinase  under  the 
control  of  albumin  promoter  (Alb-Cre)  (Supplementary  Fig.  1 1  a). 
The  tissue  specificity  of  Rpl22^^'^  incorporation  into  the  ribos¬ 
ome  was  confirmed  by  anti-HA  immunoprecipitation  of  poly¬ 
some  fractions  (Supplementary  Fig.  llb,c).  In  particular,  the 
hepatocyte-specific  transcripts  were  highly  enriched  in  the 
tagged  ribosome-associated  mRNAs  (Supplementary  Fig.  lid). 
Following  homogenization  of  liver  tissues  in  the  presence  of  LTM, 
PMY  treatment  led  to  efficient  disassociation  of  the  nonimmobi- 
lized  ribosomes  (Supplementary  Fig.  12a).  The  LTM-preserved 
ribosome  complexes  bearing  Rpl22^'^  were  enriched  by  anti-HA 
immunoprecipitation  after  RNase  I  digestion.  Metagene  analy¬ 
sis  revealed  an  elevated  ribosome  density  at  the  annotated  start 
codon,  a  result  suggesting  that  QTI-seq  worked  equally  well  in 
solid  tissue  samples  (Fig.  4a).  There  were  fewer  TISs  identified 
from  liver  cells  than  from  MEFs  (5,796  vs.  1,770),  but  they  had  a 
similar  codon  composition  ofTIS  types  (Supplementary  Fig.  13). 
In  addition,  there  was  considerable  overlap  for  the  identified  TISs 
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between  the  two  samples,  and  the  overlap  was  further  increased 
when  we  considered  only  the  aTISs  (Supplementary  Fig.  14).  The 
nonoverlapping  TISs  could  be  classified  into  differential  aTISs 
and  cell  type-specific  alternative  TISs.  For  Instance,  the  gene 
encoding  hydroxyacyl-CoA  dehydrogenase  (HADH)  showed 
a  clear  LTM  peak  at  the  aTIS  in  both  liver  and  MEF  samples 
(Supplementary  Fig.  15).  However,  an  in-frame  downstream  TIS 
was  evident  in  the  liver  sample  but  vaguely  present  in  the  MEF 
sample.  It  is  likely  that  liver  cells  can  give  rise  to  different  isoforms 
of  HADH  enzyme. 


Liver-specific  TIS  profile  in  response  to  fasting 

The  successful  application  of  QTI-seq  in  mouse  liver  offers  a 
unique  opportunity  to  explore  the  translational  response  to  a  fast¬ 
ing  condition  at  the  organismal  level.  To  this  end,  liver-specific 
RiboTag  mice  were  subject  to  overnight  fasting  followed  by  RNA- 
seq,  Ribo-seq  and  QTI-seq.  Food  deprivation  resulted  in  large  dis¬ 
assembly  of  polysomes  in  liver  tissues  (Supplementary  Fig.  12a). 
Consistently,  liver-specific  QTI-seq  uncovered  extensive  changes 
of  ribosome  density  at  the  aTISs  (Supplementary  Fig.  12b  and 
Supplementary  Table  4).  Supporting  the  quantitative  feature  of 
liver-specific  QTI-seq,  the  changes  of  aTIS  ribosome  density  in 
response  to  fasting  correlated  well  with  the  differences  of  CDS 
ribosome  occupancy  (r  =  0.512;  Fig.  4b).  Interestingly,  overnight 
fasting  elicited  distinct  translational  responses  in  liver  cells  when 
compared  to  MEFs  under  starvation  (Fig.  4c).  Although  MEFs 
showed  prominent  repression  of  ribosome  biogenesis  in  response 
to  starvation,  this  negative  translational  response  was  no  longer 
obvious  in  liver  cells  after  overnight  fasting.  In  fact,  some  ribos- 
omal  protein  genes  underwent  translational  upregulation  in 
fasted  liver.  This  observation  suggests  that,  contrary  to  the  widely 
accepted  notion,  global  inhibition  of  the  translation  machinery 
is  a  transient  response  upon  acute  nutrient  deprivation.  Under 
prolonged  fasting  conditions,  a  continuous  supply  of  ribosomal 
proteins  is  likely  needed  for  selective  protein  synthesis. 

Many  of  the  transcripts  showing  higher  aTIS  ribosome  density 
encode  components  of  the  protein  degradation  pathway,  in  par¬ 
ticular  the  proteasome  system  (Fig.  4c).  This  is  consistent  with  the 
finding  that  the  ubiquitin-proteasome  system  not  only  helps  to 
remove  misfolded  or  damaged  proteins^*  but  also  contributes  to 
intracellular  amino  acid  recycling  during  nutrient  deprivation^®. 
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Our  results — for  the  first  time,  to  our  knowledge — demonstrate 
a  translational  mechanism  that  promotes  the  proteasome  activity 
in  response  to  prolonged  fasting.  To  obtain  direct  evidence  for 
translational  control  of  the  proteasome  system,  we  cloned  the 
5'  UTR  of  several  proteasome  genes  and  fused  them  to  the  Flue 
reporter.  We  established  a  reconstituted  translation  system  pro¬ 
grammed  from  the  liver  homogenates  to  mimic  the  in  vivo  trans¬ 
lation  environment.  Using  equal  amounts  of  synthesized  mRNAs, 
we  found  that  the  presence  of  5'  UTR  of  PsmaS  or  Psmb4  resulted 
in  a  greater-than-threefold  increase  of  Flue  activity  in  liver  cells 
obtained  from  fasted  mice  but  not  in  liver  cells  from  normally  fed 
mice  (Fig.  4d).  To  ascertain  the  functional  consequence  of  this 
translational  switch  in  response  to  fasting,  we  directly  measured 
the  proteasome  activity  in  liver  homogenates  using  Proteasome- 
Glo.  Overnight  fasting  resulted  in  a  more-than-twofold  increase 
in  the  proteasome  chymotrypsin  activity  (Fig.  4e).  Our  results 
indicate  a  coordinated  regulation  between  protein  synthe¬ 
sis  and  degradation,  which  provides  an  elegant  mechanism 
for  cells  and  tissues  to  achieve  metabolic  homeostasis  under 
nutrient  deprivation. 

DISCUSSION 

Translation  initiation  is  a  crucial  point  of  regulation  in  eukaryo¬ 
tic  gene  expression,  allowing  cells  to  adapt  rapidly  to  changing 
environmental  conditions'"^.  The  significance  of  translation 
initiation  is  further  substantiated  by  the  existence  of  alternative 
translation  that  utilizes  one  or  more  potential  TISs  in  addition 
to  the  main  start  codon^^.  The  ability  to  monitor  quantitatively 
the  engagement  of  SOS  ribosomes  at  the  individual  TIS  would 
provide  an  effective  means  to  evaluate  start  codon  selection  and 
to  uncover  fundamental  principles  underlying  translational  regu¬ 
lation.  Compared  to  GTI-seq,  which  we  developed  previously^^, 
QTI-seq  represents  a  conceptually  distinct  approach  that  permits 
quantitative  profiling  of  initiating  ribosomes  in  samples  from  cells 
in  culture  and  solid  tissues. 

The  applicability  of  QTI-seq  to  solid  tissues  has  the  potential 
to  greatly  advance  our  understanding  of  the  basic  mechanisms 
of  translation  initiation  from  cellular  to  organismal  levels.  Using 
RiboTag  mice,  we  demonstrated  liver  cell-specific  profiling  of  ini¬ 
tiating  ribosomes.  This  can  serve  as  a  prototype  for  tissue-specific 
ribosome  profiling  using  other  techniques  such  as  bacTRAP^"*. 
QTI-seq  revealed  that  after  overnight  fasting,  liver  cells  had  a 
pattern  of  translational  response  distinct  from  that  of  MEFs  under 
acute  starvation,  as  follows.  First,  the  global  reduction  of  TQP 
mRNA  translation  was  no  longer  evident,  which  suggests  that  a 
continuous  supply  of  ribosomal  proteins  is  needed  for  selective 
protein  synthesis.  Second,  QTI-seq  uncovered  a  potent  transla¬ 
tional  reprogramming  for  the  proteasome  system.  The  ubiquitin- 
proteasome  system  targets  many  proteins  for  degradation,  and  the 
lethality  of  proteasome  inhibition  has  been  attributed  to  amino 
acid  scarcity^h  It  therefore  makes  intuitive  sense  that  the  upregu- 
lation  of  the  proteasome  system  is  essential  for  cell  survival  so  as 
to  continuously  supply  amino  acids.  It  will  be  of  much  interest  to 
demonstrate  this  phenomenon  in  other  types  of  tissues,  such  as 
skeletal  muscle,  in  response  to  fasting. 

METHODS 

Methods  and  any  associated  references  are  available  in  the  online 
version  of  the  paper. 


Accession  codes.  NCBI  Sequence  Read  Archive:  raw  sequence 
data  have  been  deposited  under  accession  SRA160745. 

Note:  Any  Supplementary  Information  and  Source  Datafiles  are  available  in  the 
online  version  of  the  paper. 
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Translational  control  of  the  cytosolic  stress  response  by 
mitochondrial  ribosomal  protein  LI 8 

Xingqian  Zhang',  Xiangwei  Gao*’'*,  Ryan  Alex  Coots'’^,  Crystal  S  Conn^  Botao  Liu^  &  Shu-Bing  Qian'"^ 


In  response  to  stress,  cells  attenuate  global  protein  synthesis  but  permit  efficient  translation  of  mRNAs  encoding  heat-shock 
proteins  (HSPs).  Although  decades  have  passed  since  the  first  description  of  the  heat-shock  response,  how  cells  achieve 
translational  control  of  HSP  synthesis  remains  enigmatic.  Here  we  report  an  unexpected  role  for  mitochondrial  ribosomal 
protein  LI  8  (MRPL18)  in  the  mammalian  cytosolic  stress  response.  MRPL18  bears  a  downstream  CUG  start  codon  and  generates 
a  cytosolic  isoform  in  a  stress-dependent  manner.  Cytosolic  MRPL18  incorporates  into  the  80S  ribosome  and  facilitates  ribosome 
engagement  on  mRNAs  selected  for  translation  during  stress.  MRPL18  knockdown  has  minimal  effects  on  mitochondrial  function 
but  substantially  dampens  cytosolic  HSP  expression  at  the  level  of  translation.  Our  results  uncover  a  hitherto-uncharacterized 
stress-adaptation  mechanism  in  mammalian  cells,  which  involves  formation  of  a  'hybrid'  ribosome  responsible  for  translational 
regulation  during  the  cytosolic  stress  response. 


Cell  survival  in  a  changing  environment  requires  swift  regulation  of 
gene  expression,  including  translational  control  of  existing  mRNAs'. 
Global  translation  is  generally  suppressed  in  response  to  most,  if  not 
all,  types  of  cellular  stress^’^.  However,  subsets  of  transcripts  undergo 
selective  translation  to  produce  proteins  that  are  vital  for  cell  survival 
and  stress  recovery.  One  of  the  best-known  examples  is  the  heat-shock 
protein  70  (Hsp70),  whose  synthesis  is  upregulated  in  cells  upon 
elevated  temperatures  or  exposure  to  proteotoxic  stress'*'^  Although 
transcriptional  regulation  after  heat  shock  is  well  characterized*’^, 
understanding  of  how  efficient  synthesis  of  HSPs  persists  when  the 
translation  machinery  is  generally  compromised  has  remained  elusive. 
It  is  commonly  believed  that  the  5'  untranslated  region  (5'  UTR)  of 
Hsp70  mRNA  recruits  the  translational  apparatus  in  a  cap-independent 
manner®"".  However,  neither  the  specific  translation-promoting 
features  of  Hsp70  mRNAs  nor  the  regulatory  mechanism  of  the  trans¬ 
lation  machinery  has  been  clearly  defined.  In  particular,  it  is  unknown 
whether  specialized  ribosomes  are  required  for  efficient  Hsp70 
synthesis  under  stress  conditions. 

The  ribosome  is  a  laige  ribonucleoprotein  complex  composed  of  two 
subunits  that  associate  upon  the  initiation  of  translation*^.  The  small 
subunit  decodes  mRNA,  and  the  large  subunit  catalyzes  peptide-bond 
formation.  In  mammalian  cells,  there  are  two  sets  of  ribosome  particles, 
which  reside  in  the  cytoplasm  and  the  mitochondria.  Although  ribos¬ 
omal  proteins  (RPs)  are  all  synthesized  in  the  cytoplasm,  they  assem¬ 
ble  into  functional  subunits  in  different  subcellular  compartments'®. 
The  mitochondrial  RPs  (MRPs)  are  encoded  by  nuclear  genes,  are 
synthesized  in  the  cytosol  and  are  then  imported  into  mitochondria 
for  assembly;  they  are  responsible  for  translation  of  13  mitochon¬ 
drial  niRNAs*^.  In  contrast,  cytosolic  ribosomes  are  assembled  within 


the  nucleolus  and  are  then  exported  into  the  cytoplasm  for  mRNA 
translation.  Because  of  their  distinct  composition  and  cellular  locali¬ 
zation,  there  is  believed  to  be  little  functional  connection  between 
mitochondrial  and  cytosolic  ribosomes. 

Here  we  set  out  to  investigate  whether  specialized  ribosomes  are 
required  in  the  cytosolic  stress  response.  We  report  that  MRPL18  bears 
a  hidden  CUG  start  codon  downstream  of  the  main  initiation  site. 
Stress  conditions  such  as  heat  shock  trigger  CUG-initiated  alternative 
translation,  to  generate  a  cytosolic  isoform  of  MRPL18.  We  found  that 
the  cytosolic  MRPL18  integrates  into  the  80S  ribosome  complex  in  a 
stress -dependent  manner  and  facilitates  synthesis  of  stress  proteins 
such  as  Hsp70.  Our  results  uncover  a  hitherto-uncharacterized  stress- 
adaptation  mechanism  in  mammalian  cells,  which  involves  formation 
of  a  hybrid  ribosome  that  promotes  synthesis  of  stress  proteins. 

RESULTS 

MRPL18  alternative  translation  produces  a  cytosolic  isoform 
We  used  our  previously  developed  approach,  global  translation  initia¬ 
tion  sequencing  (GTTseq),  which  allows  precise  mapping  of  alterna¬ 
tive  translation  initiation  sites  (TISs)  across  the  entire  transcriptome*®. 
The  transcript  encoding  MRPL18  bears  several  interesting  features. 
First,  MRPL18  shows  three  TISs,  with  the  annotated  start  codon  (aTIS) 
flanked  by  an  upstream  TIS  (uTIS)  and  a  downstream  TIS  (dTIS) 
(Fig.  la).  Whereas  the  uTIS  codon  is  an  AUG,  the  dTIS  uses  CUG 
as  the  initiator.  Intriguingly,  the  CUG  initiator  is  located  immedi¬ 
ately  after  the  predicted  mitochondrial  targeting  signal  (MTS)  of 
MRPL18  and  within  the  same  reading  frame.  It  is  likely  that  CUG- 
initiated  translation  produces  a  cytosolic  isoform  of  MRPL18. 
However,  the  identical  size  of  cytosolic  MRPL18  (MRPL18(cyto)) 
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Figure  1  MRPL18  undergoes  alternative 
translation  to  produce  a  cytosolic  isoform. 

(a)  Top,  schematic  of  GTI-seq  with  translation 
inhibitors  lactimidomycin  (LTM)  and 
cycloheximide  (CHX).  Raw  read  counts  per 
nucleotide  position  of  MRPL18  are  plotted 
as  a  bar  graph.  Middle,  processed  LTM  read 
density  (/?),  color-coded  by  the  corresponding 
reading  frame,  with  the  identified  TIS  positions 
marked  by  asterisks.  Bottom,  predicted  isoforms 
of  MRPL18.  (b)  Immunostaining  of  HeLa 
cells  transfected  with  MRPL18  wild  type  and 
mutants  illustrated  at  left.  Anti-myc  is  shown 
in  the  green  channel  and  MitoTracker  in  the  red 
channel.  4',6-diamidino-2-phenylindole  (DAPI) 
was  used  for  nuclear  counterstaining. 

Scale  bars,  10  pm. 


*  LTM 
mRNA 


MRPL18-myc  Anti-fnyc 


*  CHX 
mflNA  5‘ 

80  - 

>  40  - 
» 

0- 

l  20  - 

:  0  - 
X 

5  2- 
1  - 

?  0- 


o 

AAAS- 


MRPL18(WT) 


MRPLId,  NM_014161 


I j.L  ^ 
.iLliiiiuli 

■  (  il—T  '  -T~ 


MRPL18(C14A) 


Q 


MRPL18(ATG) 

ATG 


irlc 


AUG'CUG 
AUG _ 


MRPL18 

isoforms 


and  the  MTS-cleaved  mitochondrial  MRPL18 
(iMRPL18(mito))  renders  them  indistinguish¬ 
able  by  standard  immunoblotting. 

To  experimentally  confirm  the  alternative  translation  of  MRPL18, 
we  constructed  a  series  of  MRPL18  mutants  bearing  myc  tags  and 
examined  their  subcellular  localization  in  HeLa  cells.  The  wild- 
type  MRPL18  was  mostly  localized  in  mitochondria,  as  revealed  by 
immunofluorescence  staining  (Fig.  1  b).  The  vague  cytoplasmic  stain¬ 
ing  suggests  that  the  CUG-initiated  MRPL18(cyto)  was  not  a  major 
product  under  the  normal  growth  condition.  To  show  exclusively  the 
CUG-initiated  translation,  we  created  a  stop  codon  between  the  aTIS 
and  the  dTIS  to  prevent  the  synthesis  of  full-length  MRPL18.  This 
mutant,  MRPL18(C14A),  exhibited  a  clear  localization  in  the  cytosol 
and  the  nucleus  (Fig.  1  b).  We  saw  a  similar  pattern  for  MRPL18(ATG), 
a  truncated  version  of  MRPL 18  lacking  both  the  5'  UTRand  the  MTS. 
In  contrast,  MRPL18(G72A),  a  mutant  with  the  dTIS  mutated  to  CTA, 
showed  a  predominant  mitochondrial  localization.  Together,  these 
results  indicate  that  MRPL18  bears  a  hidden  downstream  CUG  start 
codon  whose  initiation  gives  rise  to  a  cytosolic  isoform  of  MRPL18. 

MRPL18  undergoes  stress-induced  alternative  translation 
Notably,  the  CUG-only  MRPL18(C14A)  had  much  lower  expression 
levels  than  other  transgenes  (Fig.  lb  and  Supplementary  Fig.  la),  a 
result  suggesting  that  the  downstream  CUG  initiator  is  less  efficient 
than  the  authentic  AUG  start  codon.  Interestingly,  MrpllS  is  a  heat 
shock-responsive  gene**,  and  a  recent  study  indicated  that  MRPL18  is 
one  of  the  direct  targets  of  heat-shock  transcription  factor  1  (HSFl)’^. 
Using  a  mouse  fibroblast  cell  line  lacking  HSFl  (ref.  18),  we  confirmed 
the  transcriptional  upregulation  of  MrpllS  upon  heat  shock 
(Supplementary  Fig.  lb).  To  investigate  whether  MRPL18  undeigoes 
translational  regulation  in  response  to  stress,  we  compared  the 
abundance  of  transfected  MRPL18  mutants  in  HeLa  cells  before  and 
after  heat  shock.  These  transgenes  maintained  similar  mRNA  levels 
in  response  to  heat  shock  (Supplementary  Fig.  Ic),  thus  allowing 
direct  evaluation  of  translational  control.  Immunoblotting  showed 
two  bands  of  transfected  MRPL18,  corresponding  to  the  MRPL18 
precursor  and  the  MRPL18  species  lacking  the  MTS  (Fig.  2a,  lane  3). 
The  transfected  MRPL18  undergoes  less  efficient  processing  than  the 
endogenous  counterpart,  presumably  owing  to  reduced  mitochondrial 
import  because  the  transgene  lacks  the  3'  untranslated  region  (3'  UTR) 
of  MRPL18  (ref.  19).  As  expected,  the  imported  exogenous  MRPL18 
incorporated  into  mitochondrial  ribosomes,  as  evidenced  by  its 
cosedimentation  with  other  mitoribosomal  proteins  (Supplementary 
Fig.  Id).  Upon  heat-shock  stress,  the  full-length  wild-type  MRPL18 
(MRPL18(WT))  showed  a  modest  increase,  whereas  we  observed 
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much  less  change  for  MRPL18(ATG),  which  lacks  both  the  5'  UTR 
and  the  MTS  region  (Fig.  2a,  lanes  9  and  10).  Remarkably,  the 
CUG-only  MRPL18(C14A)  exhibited  the  strongest  responsiveness 
despite  its  low  basal  levels  in  cells  without  stress  (Fig.  2a,  lanes  7  and  8). 
To  substantiate  this  finding  further,  we  constructed  reporters  by 
replacing  the  main  coding  region  of  MRPL18  with  firefly  luciferase 
(Flue).  Consistently  with  the  immunoblotting  results  of  MRPL18 
mutants,  the  chimeric  C14A-Fluc  e.xhibited  the  highest  increase  of 
Flue  expression  in  response  to  heat-shock  stress  (Fig.  2b). 

Having  confirmed  the  stress-induced  alternative  translation  of 
MRPL18  by  using  mutants,  we  next  sought  to  determine  whether 
MRPL18(WT)  undergoes  a  translational  switch  from  the  authentic 
AUG  to  the  downstream  CUG  in  response  to  heat-shock  stress. 
It  is  challenging  to  unequivocally  monitor  the  newly  synthesized 
MRPL18(cyto)  because  of  the  high  basal  levels  of  MRPL18(mito) 
present  in  cells  before  stress.  To  circumvent  this  limitation,  we  engi¬ 
neered  a  reversible  destabilization  domain  (DD)  by  fusing  it  to  the 
COOH  terminus  of  MRPL18(WT)  (Fig.  2c).  This  system  prermits  tem¬ 
poral  examination  of  newly  synthesized  MRPL18-DD  after  addition  of 
Shield- 1,  a  cell-permeable  drug  that  binds  the  DD  to  protect  the  protein 
from  degradation^®.  HeLa  cells  transfected  with  plasmids  expressing 
MRPL18-DD  e.xhibited  minimal  anti-myc  signals  in  the  absence  of 
Shield-1  (Supplementary  Fig.  2a).  Treatment  with  Shield-1  stabilized 
MRPL18-DD,  which  had  a  clear  mitochondrial  localization  (Fig.  2c  and 
Supplementary  Fig.  2b).  Remarkably,  heat-shock  stress  before  Shield- 1 
addition  resulted  in  substantial  anti-myc  signals  in  the  cytosol  as  well  as 
the  nucleus  (Fig.  2c  and  Supplementary  Fig.  2c).  This  observation  was 
not  due  to  mitochondrion  leakage,  because  the  stress- induced  Hsp60 
remained  exclusively  in  mitochondria  after  heat  shock  This  result  con¬ 
firms  that  heat-shock  stress  triggers  CUG- mediated  alternative  transla¬ 
tion  within  the  wild-type  sequence  context  of  MRPL18. 

MRPL18(cyto)  is  dependent  on  elF2a  phosphorylation 
What  is  the  mechanism  underlying  stress- induced  alternative  trans¬ 
lation  of  MRPL18?  We  looked  into  eukaryotic  initiation  factor  2a 
(eIF2a),  whose  phosphorylation  regulates  ATF4,  in  a  classical  example 
of  alternative  translation  triggered  by  many  stress  conditions^'"^^. 
Consistently  with  a  previous  report^^,  heat-shock  stress  also  triggered 
eIF2a  phosphorylation  (Suppilementary  Fig.  3a).  To  address  the  role 
of  eIF2a  phosphorylation  in  the  alternative  translation  of  MRPL18, 
we  transfected  MRPL18-Fluc  reporters  into  a  mouse  embryonic  fibro¬ 
blast  (MEF)  cell  line  bearing  a  nonphosphorylatable  eIF2a  in  which 


NATURE  STRUCTURAL  &  MOLECULAR  BIOLOGY  VOLUME  22  NUMBERS  MAY  2015 


405 


ARTICLES 


MRPL18 


Plasmid  X  ^  ^  xC> 

<?  ^  MW 

HS-  +  -  +  -  +-  +-  +  (kDa) 

cc 

2 


|1 
OC  *“ 


P 

25 

09 

1  2  3  4  5  6  7  6 

’  □  No  HS 

-  ■  Post  HS  1 

ni  nl  nl 

9  10 

MRPL18(WT)  OD 


MRPL18  (mtto)  (cyto) 

^  <3*^ 

X^-ShleW-1 
t  ‘ 

CTG  Myc 

MRPL18  (mito)  (cyto) 
Anti-Hsp60  Anti-myc  '-At- 


Post  HS 
K  Shdlid-1 


MRPt.18-Fluc 

CTG 

CTA 

V 

I  rrrft 


□  No  HS 
■  Post  HS 


r^i 


ATG 


□  NoHS 
El  Post  HS 


\.+St>ieW-1 

-<M> 


MRPt18>RiiC  elF2a(S/S) 


elF2a  (A/A) 


— I — I — I — I — 

0  1  2  3  4  5  (lotd) 


Figure  2  MRPL18  undergoes  alternative 
translation  in  a  stress-dependent  manner. 

(a)  Representative  immunoblot  results  of  HeLa 
cells  transfected  with  the  indicated  plasmids 
and  collected  before  heat  shock  (HS)  or  2.5  h 
after  heat  shock  (43  X.  1  h).  Top  bands  are 
uncleaved  MRPL18  precursors.  (J-actin  is  a 
loading  control.  The  bottom  graph  shows  the 
relative  levels  of  transgenes  (small  isoform) 
normalized  to  p-actin.  Error  bars,  s.e.m. 

(n  =  3  cell  cultures,  *P<  0.01  by  two-tailed 
Student's  ftest).  MW,  molecular  weight; 
exp,  exposure,  (b)  Fiuc  reporter  assays  in  HeLa 
cells  transfected  with  MRPL18-Fluc  constructs 
shown  at  left.  Flue  activities  after  heat  shock 
are  normalized  to  those  in  nonstressed  cells. 
Error  bars,  s.e.m.  (n  =  3  cell  cultures). 

(c)  Immunostaining  of  HeLa  cells  transfected 
with  MRPL18-DD  with  or  without  heat  shock 
in  the  presence  of  Shield-1.  Anti-Hsp60 
antibody  is  shown  in  the  red  channei  and 
anti-myc  in  the  green  channel.  Scale  bars, 

10  pm.  A  schematic  of  MRPL18-DD  fusion 
protein  in  the  absence  or  presence  of  Shield-1 
is  shown  at  top.  (d)  Flue  reporter  assays  in 
elF2a  (S/S)  and  elF2a  (/V/A)  MEFs  transfected 
with  MRPLIS-Fluc  constructs  shown  in  b. 

Flue  activities  after  heat  shock  are  normalized 
to  those  of  nonstressed  cells.  Error  bars, 
s.e.m.  (n  =  3  cell  cultures).  Original  blot  images 
are  in  Supplementary  Data  Set  1. 


SerSl  (S/S  genotype)  was  mutated  to  an  alanine  (A/A  genotype)^*. 
As  expected,  wild-type  MEF  (S/S)  cells  showed  a  similar  pattern  of 
Flue  expression  as  HeLa  cells  (Fig.  2d).  In  particular,  the  CUG-only 
MRPL18(C14A)  exhibited  the  highest  increase  of  Flue  (more  than 
four-fold)  in  response  to  heat-shock  stress,  whereas  the  expression 
level  of  the  AUG-only  MRPL18(G72A)  remained  unchanged. 
Remarkably,  the  stress  responsiveness  of  MRPL18(C14A)  was  com¬ 
pletely  abolished  in  MEF  (A/A)  cells  (Fig.  2d  and  Supplementary 
Fig.  3b),  thus  indicating  that  the  alternative  translation  of  MRPL18 
is  dependent  on  eIF2a  phosphorylation. 

MRPL18(cyto)  integrates  into  cytosolic  SOS  ribosomes 

Given  that  the  authentic  role  of  MRPL18  is  to  constitute  the  ribosome 

complex  within  mitochondria,  we  speculated  that  MRPL18(cyto)  may 


integrate  into  the  cytosolic  80S  ribosome  complex  under  stress  condi¬ 
tions.  To  test  this  hypothesis,  we  assessed  the  behavior  of  endogenous 
MRPL18  as  well  as  transfected  mutants  in  HeLa  cells.  We  first  purified 
cytosolic  ribosome  complexes  with  affinity  immunoprecipitation  (IP) 
in  order  to  eliminate  mitochondrial  ribosome  contamination.  We  also 
converted  the  polysomes  into  monosomes  by  RNase  I  digestion  before 
IP  to  exclude  indirect  pulldown  of  RNA-binding  proteins  (Fig.  3a). 
Endogenous  MRPL18,  but  not  other  mitochondrial  proteins,  was  read¬ 
ily  precipitated  from  stressed  cells  by  either  anti-RPL4  or  anti-RPS6 
antibodies  (Fig.  3a,  lanes  2  and  6).  In  addition,  we  recovered  endog¬ 
enous  MRPLI8  from  the  polysome  fraction  after  heat-shock  stress 
and  found  that  this  redistribution  was  sensitive  to  the  translation 
inhibitor  cycloheximide  (CHX)  (Supplementary  Fig.  4).  Therefore, 
the  cytosolic  ribosome-associated  MRPL18  is  newly  synthesized  after 
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Figure  3  Cytosolic  MRPL18  incorporates  into  the 
80$  ribosome  in  a  stress-dependent  manner. 

(a)  Detection  of  MRPL18  in  endogenous 
ribosomes  purified  from  HeLa  cells  with  or 
without  heat  shock.  Immunoblotting  of  RNase 
l-digested  cell  lysates  immunoprecipitated  with 
anti-RPL4  (left)  or  anti-RPS6  (right)  antibodies. 

Input  and  immunoprecipitates  (IP)  are  indicated. 

Cyto,  cytoplasmic;  mito,  mitochondrial. 

(b)  Immunoblot  detection  of  myc-tagged 
MRPL18(ATG)  from  endogenous  ribosomes 
isolated  by  sucrose  cushion  from  HeLa  cells  with 
or  without  heat  shock.  HeLa  cells  transfected 
with  GFP  are  a  control.  Total  lysate  and  ribosome 
(ribo)  pellet  are  indicated,  (c)  Immunoblot 
detection  of  cytosolic  ribosomal  proteins  from 
anti-myc  immunoprecipitates.  Samples  are  HeLa 
cells  transfected  with  myc-tagged  MRPL18(ATG), 

with  or  without  heat  shock,  and  HeLa  cells  transfected  with  GFP  and  MRPL18(G72A)  as  controls.  Cell  lysates  were  treated  with  RNase  I  before  anti-myc  IP. 
Throughout  figure,  p-actin  is  a  loading  control.  Experiment  schematics  are  shown  at  top.  Original  blot  images  are  in  Supplementary  Data  Set  1. 
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Figure  4  Cytosolic  MRPL18  incorporates  into  the  SOS  ribosome  in  a 
phosphorylation-dependent  manner,  (a)  Detection  of  phosphorylated 
MRPL18  in  HeLa  cells  with  or  without  heat-shock  stress.  Immunoblotting 
of  whole  cell  lysates  with  or  without  A.-phosphatase  (X-PP)  treatment, 
separated  with  Phos-tag  acrylamide  gel  (MRPLIS  only)  or  standard 
PAGE.  CHX  treatment  during  heat  shock  was  included  to  demonstrate 
phosphorylation  of  newly  synthesized  MRPLIS.  (b)  Detection  of 
phosphorylated  (P)  MRPLIS  in  HeLa  cell  fractions  with  or  without  heat- 
shock  stress.  Immunoblotting  of  HeLa  lysates  with  or  without  heat  shock, 
fractionated  with  a  mitochondria  fractionation  kit  and  separated  with 
Phos-tag  acrylamide  gel  (MRPLIS  only)  or  standard  PAGE,  (c)  Detection 
of  phosphorylated  MRPLIS  in  ribosome  fractions  of  HeLa  cells  with  or 
without  heat-shock  stress.  Immunoblotting  of  total  lysate  and  ribosome 
fractions  prepared  by  sucrose  cushion,  separated  with  Phos-tag  acrylamide 
gel  (MRPLIS  only)  or  standard  PAGE,  (d)  Detection  of  phosphorylated 
MRPLIS  in  ribosomes  purified  from  HeLa  cells  with  or  without  heat-shock 
stress.  Immunoblotting  of  RNase  l-digested  total  lysates  and  anti-RPL4 
immunoprecipitates,  separated  with  Phos-tag  acrylamide  gel  (MRPLIS 
only)  or  standard  PAGE.  Throughout  figure,  P-actin  is  a  loading  control. 
Original  blot  images  are  in  Supplementary  Data  Set  1. 


stress.  We  further  confirmed  the  incorporation  of  MRPL18(cyto) 
into  the  SOS  ribosome  complex  by  using  transfected  MRPL18(ATG). 
Both  cosedimentation  binding  analysis  and  anti-myc  IP  revealed 
a  robust  association  of  MRPL18(ATG)  with  cytosolic  ribosomes 
in  a  stress-dependent  manner  (Fig.  3b,  lane  8  and  Fig.  3c,  lane  6). 
In  contrast,  the  AUG-only  MRPLl  8(G72A)  did  not  coprecipitate  with 
any  cytosolic  RPs  before  or  after  heat-shock  stress  (Fig.  3c). 
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Phosphorylated  MRPLl  8(cyto)  integrates  into  BOS  ribosomes 
Interestingly,  very  little  transfected  MRPL18(ATG)  was  associated 
with  the  805  ribosome  in  unstressed  cells  in  spite  of  comparable 
expression  levels.  To  investigate  the  mechanism  of  stress- induced 
incorporation,  we  examined  the  phosphorylation  status  of  MRPLIS 
by  using  a  Phos-tag  acrylamide  gel.  An  MRPLIS  band  with  slower 
migration  was  clearly  discernible  upon  heat-shock  stress,  and  it  was 
sensitive  to  phosphatase  treatment  (Fig.  4a).  Addition  of  the  transla¬ 
tion  inhibitor  CHX  completely  abolished  this  species,  thus  suggest¬ 
ing  that  newly  synthesized  MRPLl  8  is  phosphorylated  during  stress. 
Importantly,  we  found  the  phosphorylated  MRPLl  8  exclusively  in  the 
cytosolic  fraction  (Fig.  4b),  and  we  recovered  only  the  phosphorylated 
species  from  the  cytosolic  80S  ribosome  (Fig.  4c,d).  We  next  searched 
for  the  kinase  responsible  for  MRPLl  8  phosphorylation.  Interestingly, 
MRPL18  was  predicted  to  be  a  substrate  of  the  tyrosine  protein  kinase 
Lyn^^.  Indeed,  both  the  Lyn-specific  chemical  inhibitor  bafetinib 
(Fig.  5a)  and  small  hairpin  RNA  (shRNA)-mediated  knockdown 
of  Lyn  (Fig.  5b)  reduced  the  phosphorylation  level  of  MRPLIS. 
Notably,  Lyn  is  a  member  of  the  Src  family,  whose  kinase  activity  is 
increased  upon  heat-shock  stress^*.  These  results  collectively  indicate 
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that  both  the  production  and  function  of  MRPL18(cyto)  are  under 
tight  control  in  response  to  stress. 

MRPLl  8(cyto)  promotes  stress-protein  synthesis  during  stress 
The  stress-inducible  feature  of  MRPL18( cyto)  is  suggestive  of  its  regu¬ 
latory  role  in  the  cytosolic  stress  response.  To  elucidate  its  physiological 
function,  we  knocked  down  MRPL18  in  HeLa  cells  with  SMARTpool 
small  interfering  RNAs  (siRNAs).  Reduction  in  MRPL18  expression 
had  minimal  effects  on  cell  growth  and  did  not  alter  the  rate  of  global 
protein  synthesis  (Supplementary  Fig.  5a).  In  addition,  we  observed 
minimal  effects  of  MRPL18  depletion  on  mitochondrial  transla¬ 
tion  after  short-term  siRNA  knockdown  or  long-term  lentivirus- 
based  shRNA  knockdown  (Supplementary  Fig.  5b,c).  However, 
upon  heat-shock  treatment,  the  induction  of  Hsp70  expression  in  the 
cytosol  was  largely  dampened  in  cells  with  MRPL 1 8  knockdown  com¬ 
pared  to  cells  transfected  with  control  siRNA  (Fig.  6a).  To  substanti¬ 
ate  this  finding  further,  we  applied  shRNA  to  HeLa  cells  and  MEFs 
(Fig.  6b  and  Supplementary  Fig.  6).  shRNA-mediated  MRPL18 
knockdown  resulted  in  impaired  Hsp70  induction  in  both  cells  after 
heat-shock  stress.  This  result  was  not  due  to  transcriptional  deficiency 
of  Hsp70  gene  expression.  In  fact,  MEFs  with  MRPLI8  knockdown 
demonstrated  even  higher  Hsp70  transcript  levels  than  the  control 
cells  in  response  to  heat-shock  stress  (Supplementary  Fig.  6d). 


Figure  5  Cytosolic  MRPLIS  undergoes  Lyn-nnediated  phosphorylation. 

(a)  Detection  of  phosphorylated  MRPLIS  in  HeLa  ceils  in  the  presence 
of  increasing  doses  of  bafetinib  immediately  after  heat-shock  stress. 
Immunoblotting  of  whole  cell  lysates  separated  with  Phos-tag  acrylamide 
gel  (MRPLIS  only)  or  standard  PAGE,  (b)  Detection  of  phosphorylated 
MRPLIS  in  HeLa  cells  transfected  with  SMARTpool  siRNAs  targeting 
Lyn  or  control  shRNA  targeting  scrambled  sequences.  Immunoblotting  of 
whole  cell  lysates  separated  with  Phos-tag  acrylamide  gel  (MRPLIS  only) 
or  standard  PAGE.  Throughout  figure,  ft-actin  Is  a  loading  control.  Original 
blot  images  are  in  Supplementary  Data  Set  1. 
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Figures  Cytosolic  MR  PL  18  promotes 
Hsp70  biosynthesis  after  heat-shock  stress. 

(a)  Examination  of  heat  shock-induced  Hsp70 
synthesis  in  HeLa  cells  transfected  with 
SMARTpool  siRNAs  targeting  MRPL18  or  control 
shRNA  targeting  scrambled  sequences. 

Immunoblottingof  cell  lysates  collected  at 
indicated  times  after  heat  shock.  N,  no  heat  shock. 

(b)  Examination  of  heat  shock-induced  Hsp70 
synthesis  in  MEFs  infected  with  lentiviruses 
expressing  shRNAs  targeting  MRPL18  or 
control  shRNA.  Immunoblotting  of  cell  lysates 
collected  at  the  indicated  times  after  heat 
shock,  (c)  Examination  of  heat  shock-induced 
Hsp70  synthesis  in  elF2a  (S/S)  and  elF2a  (A/A) 

MEFs  after  heat-shock  stress,  (d)  Examination 
of  polysome-enriched  Hsp70  mRNA  in  MEFs 
infected  with  lentiviruses  expressing  shRNA  targeting  MRPL18  or  control  shRNA.  Quantitative  PCR  measuring  Hsp70  and  p-actin  mRNA  levels  in  total 
RNA  extracted  from  polysome  fractions.  Relative  levels  are  normalized  to  the  total.  OD254,  optical  density  at  254  nm.  Error  bars,  s.e.m.  (n=  3  cell  cultures, 
*P<  0.05  by  two-tailed  Student's  f  test).  Throughout  figure,  P-actin  is  a  loading  control.  Original  blot  images  are  in  Supplementary  Data  Set  1. 
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Given  that  the  de  novo  synthesis  of  MRPL18(cyto)  relies  on  alterna¬ 
tive  translation,  we  reasoned  that  suppressing  CUG  initiation  would 
have  a  similar  effect  as  MRPL18  knockdown.  Indeed,  stress- induced 
Hsp70  synthesis  was  substantially  reduced  in  eIF2a  (A/ A)  cells,  whose 
eIF2a  is  not  phosphorylatable  (Fig.  6c).  Because  the  incorporation  of 
MRPL18(cyto)  into  the  80S  ribosome  is  phosphorylation  dependent, 
we  examined  the  effects  of  reduced  Lyn  kinase  activities  in  the  stress 
response.  Treatment  with  the  Lyn  inhibitor  bafetinib  or  Lyn  knock¬ 
down  dramatically  suppressed  Hsp70  induction  after  heat-shock 
stress  (Fig.  5a,b).  The  deficient  Hsp70  synthesis  is  specific  to  MRPL18 
because  knocking  down  other  mitochondrial  ribosomal  proteins  did 
not  influence  the  cytosolic  stress  response  at  all  (Supplementary 
7a,b).  These  data  collectively  highlight  MRPL18(cyto)  as  a 
critical  mediator  of  the  cytosolic  stress  response. 


MRPL18(cyto)  promotes  ribosome  engagement  on  stress  mRNAs 
The  participation  of  MRPL18(cyto)  in  Hsp70  synthesis  suggests  that 
a  specialized  ribosome  might  be  required  for  mRNA  translation 
under  stress  conditions.  MRPLl  8  is  a  homolog  of  cytoplasmic  RPL5, 
according  to  sequence  alignment.  Much  like  RPL5,  MRPL18  binds  to 
5S  rRNA  and  is  believed  to  import  the  5S  rRNA  into  mitochondria^*. 
An  in  vitro  RNA  binding  assay  confirmed  that  5S  rRNA  strongly  asso¬ 
ciates  with  recombinant  MRPL18(ATG)  but  not  GFP  (Supplementary 
Fig.  7d,e).  However,  we  found  little  evidence  indicating  that  MRPL18 
replaces  RPL5  in  binding  to  the  80S  ribosome  (Fig.  3).  In  fact,  RPL5 
knockdown  also  reduced  Hsp70  synthesis  after  heat-shock  stress 
(Supplementary  Fig.  7c).  Notably,  MRPL18  does  not  seem  to  interact 
with  separated  40S  or  60S  subunits  (Supplementary  Fig.  7f).  It  is 
thus  likely  that  MRPL18(cyto)  serves  as  an  extra  RP  in  binding  to  the 


Figure  7  Cytosolic  MRPL18  is  essential 
for  induced  thermal  tolerance,  (a)  Thermal 
tolerance  analysis  of  MEFs  lacking  MRPL18. 
MEFs  transfected  with  shRNA  targeting 
MRPL18  or  control  shRNA  targeting  scrambled 
sequences  were  primed  with  mild  heat  shock 
(43  °C  for  30  min)  and  allowed  to  recover  at 
37  “C  for  5  h  (pre).  Severe  heat  shock  (45  °C 
for  1  h;  lethal)  was  then  applied  before  the  cell 
viability  assay.  Error  bars,  s.e.m.  (n  =  4  cell 
cultures,  *P<  0.01  by  two-tailed  Student’s 
f  test).  Control,  no  heat  shock,  (b)  Examination 
of  apoptotic  markers  in  MEF  cells  as  in  a  by 
immunoblotting  with  the  indicated  antibodies. 

(c)  Thermal  tolerance  assays  of  MEFs  after 
reintroduction  of  Fisp70.  MEFs  as  in  a  were 
infected  with  recombinant  adenoviruses 
(AdV)  encoding  GFP  or  Hsp70  before 
thermal  tolerance  analysis.  Error  bars,  s.e.m. 
{n=  3  cell  cultures,  *P<  0.01  by  two-tailed 
Student’s  f  test).  Control,  uninfected  cells. 

(d)  A  model  for  MRPL18  in  translational  control 
during  stress  conditions.  Stress  conditions 
such  as  heat  shock  trigger  trimerization  and 
activation  of  HSFl  (1),  which  turns  on  genes 
including  MrpIlS and  Hspala.  In  nonstressed 
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cells,  MRPL18  translation  is  initiated  from  the  annotated  AUG,  generating  MRPL18  containing  a  mitochondria  localization  signal  (blue).  In  stressed 
cells,  MRPL18  undergoes  alternative  translation  from  the  downstream  CUG  start  codon,  which  results  in  a  cytoplasmic  isoform  of  MRPL18  (red)  (2). 
MRPL18(cyto)  incorporates  into  the  805  ribosome  complex,  facilitating  the  engagement  of  mRNAs  highly  expressed  under  stress,  such  as  Hspala 
(green)  (3).  The  efficient  synthesis  of  heat-shock  proteins  contributes  to  stress  recovery.  Original  blot  images  are  in  Supplementary  Data  Set  1. 
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80S  ribosome.  5S  rRNA  forms  the  central  protuberance  of  the  lai^e 
ribosomal  subunit^®.  This  region  is  of  particular  importance  because 
it  nears  the  place  where  the  ribosomal  large  and  small  subunits, 
the  decoded  mRNA  and  the  peptidyl  tRNA  all  come  together.  We  pos¬ 
tulate  that  the  presence  of  an  extra  RP  promotes  SOS  ribosome  engage¬ 
ment  on  mRNAs  during  stress  conditions.  To  test  this  hypothesis,  we 
examined  the  distribution  of  Hsp70  mRNA  in  the  polysome  fractions 
of  stressed  MEF  cells  with  or  without  MRPL18  knockdown  (Fig.  6d). 
With  comparable  total  mRNA  levels,  the  FIsp70  transcript  showed 
much  less  enrichment  in  the  polysomes  of  MEFs  lacking  MRPL18. 
In  contrast,  control  |i-actin  mRNA  exhibited  only  minor  reduction 
in  the  polysomes  of  these  cells,  presumably  because  of  the  delayed 
stress  recovery  as  a  result  of  impaired  Hsp70  synthesis.  Indeed, 
MEFs  lacking  MRPL18  showed  less  polysome  formation  during  stress 
recovery  (Fig.  6d). 

MRPL18  is  essential  for  induced  thermotolerance 
Because  of  the  essential  role  of  chaperone  molecules  in  cell  survival, 
we  predict  that  cells  lacking  MRPLl  8  should  have  attenuated  thermo¬ 
tolerance.  Induced  thermotolerance  allows  cells  to  survive  a  normally 
lethal  temperature  if  they  are  first  conditioned  at  a  milder  temperature. 
Without  preconditionir^,  MEFs  with  or  without  MRPLl  8  knockdown 
were  equally  susceptible  to  severe  heat  stress  at  45  °C  for  1  h  (Fig.  7a). 
Preexposure  at  43  °C  for  30  min  resulted  in  a  nearly  complete  pro¬ 
tection  of  control  MEF  cells  from  severe  heat  stress.  However,  we  no 
longer  observed  the  induced  thermotolerance  in  MEFs  with  MRPLl  8 
knockdown.  As  expected,  the  reduced  cell  viability  in  the  absence  of 
MRPL18  was  laigely  due  to  increased  apoptosis  (Fig.  7b).  In  agree¬ 
ment  with  the  suppressed  chaperone  biosynthesis  in  cells  lacking 
MRPL18,  reintroduction  of  Hsp70  by  recombinant  adenoviruses 
completely  restored  the  thermotolerance  (Fig.  7c).  Given  the  broad 
range  of  chaperone  function  in  cell  physiology,  it  would  be  of  much 
interest  to  investigate  possible  protective  roles  of  MRPLl 8  against 
cellular  stressors  beyond  heat  shock. 

DISCUSSION 

Cells  in  nearly  all  living  organisms  respond  to  heat-shock  stress  by 
marked  transcriptional  alterations  and  rapid  translational  repro¬ 
gramming’.  Despite  severe  inhibition  of  the  translation  machinery 
under  stress  conditions,  efficient  synthesis  of  stress  proteins  persists 
by  a  mechanism  that  has  not  been  completely  understood.  Here,  we 
uncovered  a  molecular  mechanism  underlying  the  active  translation 
of  mRNAs  highly  expressed  during  stress.  We  show  that  heat-shock 
stress  triggers  a  previously  unrecognized  alternative  translation  of 
MRPLl  8,  which  generates  a  cytoplasmic  isoform  of  the  mitochondrial 
ribosomal  protein.  Remarkably,  the  cytoplasmic  version  of  MRPL18 
incorporates  into  the  SOS  ribosome  complex  in  a  stress-dependent 
manner.  The  stress-induced  formation  of  specialized  ribosomes, 
together  with  the  cap- independent  translation  initiation  mechanism, 
ensures  efficient  translation  of  mRNAs  under  nonfavorable  conditions 
(Fig.  7d).  Our  results  thus  provide  a  new  paradigm  for  translational 
regulation  under  stress  conditions,  which  involves  ribosome  speciali¬ 
zation  after  an  initial  alternative  translation  event. 

A  growing  body  of  evidence  has  suggested  that  ribosome  hetero¬ 
geneity  prevails  across  species,  under  different  developmental  stages 
and  in  various  tissues^’’^^.  Variation  in  ribosome  composition,  in  both 
rRNA  and  RPs,  provides  a  regulatory  mechanism  to  the  translation 
machinery.  A  clear  example  is  illustrated  in  Escherichia  coli,  in  which 
the  stress-induced  endonuclease  MazF  cleaves  the  16S  rRNA  and 
removes  the  anti-Shine-Dalgarno  sequence^^.  The  resultant  ‘stress 
ribosome’  selectively  translates  the  leaderless  mRNAs,  a  group  of 


transcripts  also  generated  by  MazF.  In  eukaryotes,  certain  RPs  have 
been  found  to  control  transcript  selectivity  during  translation.  For 
example,  RPL38  is  required  for  translation  of  homeobox  mRNAs^'*, 
whereas  RPL40  appears  to  control  translation  of  vesicular  stomatitis 
virus  mRNAs^^.  In  this  report,  we  discovered  in  mammalian  cells  an 
unusual  mechanism  that  allows  HSP  mRNA  translation  to  escape 
from  the  shutoff  of  global  protein  synthesis.  Instead  of  altering  the 
rRNA  structure,  the  mammalian  stress  ribosome  uses  an  extra  RP  that 
functions  in  mitochondria  at  basal  state.  MRPL18(cyto)  incorporation 
may  alter  ribosome  conformation  and/or  stabilize  ribosome  engage¬ 
ment  on  mRNAs  that  are  highly  expressed  under  stress  conditions. 
Alternatively,  the  presence  of  MRPL18(cyto)  may  cause  recruitment 
of  additional  factors  facilitating  initiation  and  elongation.  In  any  case, 
stress-induced  ribosome  heterogeneity  permits  translational  repro¬ 
gramming  without  rebuilding  the  entire  translational  machinery. 

Another  interesting  phenomenon  revealed  by  our  data  is  the  func¬ 
tional  connection  between  mitochondrial  and  cytoplasmic  ribosomes. 
The  55S  mitochondrial  ribosome  differs  substantially  from  the 
80S  ribosome  in  eukaryotic  cytoplasm  and  the  70S  ribosome  in 
prokaryotes.  In  mammalian  cells,  the  mitochondrial  ribosome  complex 
has  a  higher  content  of  protein  than  rRNA  components*'*'^.  Althoi^h 
many  MRPs  are  distinctive  and  evolving  rapidly,  MRPL18  has  a  close 
homolog  in  E.  coli.  Unlike  many  other  RPs,  MRPL18  has  apparently 
evolved  to  become  a  stress-inducible  gene  in  mammals.  The  posses¬ 
sion  of  an  alternative  translation  feature  further  renders  MRPL18  a 
critical  regulator  of  the  stress  response.  We  conclude  that  MRPL18 
actively  participates  in  stress  adaptation  by  potentiating  the  cellular 
translation  machinery  to  achieve  a  robust  cytosolic  stress  response. 

METHODS 

Methods  and  any  associated  references  are  available  in  the  online 
version  of  the  paper. 

Note:  Any  Supplementary  Information  and  Source  Data  files  are  available  in  the  online 
version  of  the  paper. 
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